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ABSTRACT 

Context. Young stars are formed within dusty discs. The grains in the disc are originally of the same size as interstellar dust, around 
0.1 yum. Models predict that these grains will grow in size through coagulation. Observations of the silicate features at fjm wavelengths 
are consistent with growth from submicron to micron sizes in selected sources whereas the slope of the SED at longer wavelengths 
traces growth up to mm sizes and larger. 

Aims. We here look for a correlation between these two grain growth indicators. 

Methods. A large sample of T-Tauri and Herbig-Ae/Be stars, spread over five different star-forming regions, was observed with the 
Spitzer Space Telescope at 5-13 yum; a subsample was observed at mm wavelengths. We complement this subsample with data from 
the literature to maximise the overlap between lum and mm observations and search for correlations. Synthetic spectra are produced 
to determine which processes may produce the dust evolution. 

Results. Dust disc masses in the range < 1 to 7 x 10"'* Mq are obtained. Most sources have a mm spectral slope consistent with 
grain growth. There is a tentative correlation between the 10-yum silicate feature and the slope of the SED between 1 and 3 mm. The 
observed sources seem to be grouped per star-forming region in the yum-vs-mm diagram. The modelling results show that first the 
10-yum feature becomes flatter and subsequently the mm slope becomes shallower. Grain size distributions shallower than that of the 
ISM and/or bright central stars are required to explain specific features. Settling of larger grains towards the disc midplane affects 
10-yum feature, but hardly the mm slope. 

Conclusions. The tentative correlation between the strength of the 10-//m feature and the mm slope suggests that the inner and outer 
disc evolve simultaneously. Dust with a mass dominated by ~mm-sized grains is required to explain the shallowest mm slopes. Other 
processes besides grain growth may also be responsible for the removal of small grains. Observations with future, more sensitive 
telescopes are required to provide the necessary statistics to study these processes of disc and dust evolution. 

Ksy words, circumstellar matter - planetary systems: protoplanetary discs - stars: formation 



1. Introduction 

A long-standing problem in planet formation is how tiny inter- 
stellar dust particles of less than a micron in size coagulate and 
grow to eventually form planets, thousands of kilometres in size. 
It is in the very nature of this field that it has to be studied at 
various levels, since different physical processes dominate dur- 
ing the various phases. The first steps, which lead to dust grains 
of about a decimetre in size, are studied both in the laboratory 
and with computer simulations (see Domi nik et al.|2007 , Blum. 
|& Wurm]|2008[ for detailed reviews). Local concentrations of 
boulders and subsequent gravitational collapse may then lead to 
the formation of planetesimals several hundreds of kilometres in 
size (e.g., |Johansen et al.|2007[ ). This paper focuses on the obser- 
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vational signatures of (sub)micron-sized grains up to centimetre- 
sized pebbles. 

The InfraRed Spectrograph (IRS) on-board the Spitzer Space 
Telescope has provided a wealth of mid-infrared (5^0 fim) 
spectra from discs around pre-main-sequence stars (e.g. |Kessler- 
|Silacci et al.|2006 Furlan et al.|2006^ . The spectra of these ob- 
jects are often dominated by silicate emission features at 10 and 
20 fim. In young stellar objects, these features are formed in the 
upper atmosphere of the hot inner disc. The varying strength and 
the shape of these features can be naturally explained by differ- 
ent grain sizes in the upper layers of the inner disc, with strong, 
pointed features being representative of ~0. 1 /vm-sized grains 
and flatter features coming from dust grains of several yum in 
size (Kessler-Silacci et al.|[2006 ). These results confirm earlier 
results from the Infrared Space Observatory (Bouwman et al. 
|2001[ |Meeus et al.pOOl) and from ground-based observations 
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(e.g., |Przygodda et al.|2003| l. It has been suggested that crystalU- 
sation has a similar effect on the 10-jum feature as grain growth 
(e.g., Honda et al. 2003, Meeus et al. 2003 ). However, this effect 
is minimal and the dominating factor for the strength and shape 



of the IQ-fim feature is the characteristic grain size (Olofsson 
let al.|2009| l. 

Because the 10-/vm feature only probes the surface layers 
of the inner disc, a stronger, more peaked feature could also 
be due to the settling of larger, micron-sized grains towards the 
mid-plane. As the larger grains settle and the small ones remain 
suspended in the upper layers, the surface becomes dominated 
by small grains, creating a strong silicate band. |Dullemond &| 
[Dominik^ (2008^ 1 investigated this idea through theoretical mod- 
els. They find that settling can in principle explain the different 
shapes of the 10-jum feature, but only in quite specific cases, 
so that overall grain growth is still the most likely explanation 
for the flattening of these features. Recent interferometric ob- 
servations of the IQ-fim spectral region in discs around ~1 and 
2-3-Mo objects show that the grains closer to the central star 
are both larger and more crystalline than those further out in the 
disc (see, e.g., the recent review by , van Boekel 2008 ). Hence, the 
evolution of the lO-fim feature may be caused by a combination 
of grain growth and crystallisation and appears to progress from 
the inner disc outwards. On the other hand, analysis of the longer 
wavelenth mid-infrared crystalline features indicates significant 
growth and crystallisation in the outer disc as well (Olofsson, 
let al.|2009| l. 

Whereas the mid-infrared region potentially provides infor- 
mation on the growth of grains from interstellar, submicron sizes 
to sizes of several microns, the growth to larger sizes can only 
be probed by submillimetre (submm), millimetre (mm), and oc- 
casionally centimetre (cm) observations. Ground-breaking work 
was done by Beckwith et al.| ( T990[ l and Beckwith & Sargent 
( 1991|l, both analytically studying the emission of dust grains 



and obtaining the first submm slopes by observing a large sam- 
ple of young stellar objects at mm wavelengths. More recently, 
Andrews & Wilhams ( 2005[ l performed a sensitive single-dish 



submm continuum survey of 153 young stellar objects in the 
Taurus-Auriga star-formation region, including a large amount 
of archival and literature data. They found that the submm 
slope between 350 fj.m and 1.3 mm could be well described by 
a = 2.0 ± 0.5, where Fy oc y", while the value for the interstel- 
lar medium is a ^ 3.7 (cf. l3raine'2 006| l. [Andrews & Williams| 
( |200 5 ) interpreted this shallow slope as a combined effect of a 
contribution from optically thick regions in the disc and grain 
growth. It should be noted, however, that the sources in this study 
were spatially unresolved, and the (sub)mm emission may have 
a significant contribution from surrounding (envelope) material. 
More recently, interferometric studies of several dozen T-Tauri 
stars gave values of a < 3.0 ([Rodmann et al.||2006| | Andrews] 
|& Williams||2007| |Lommen et al.||2007| ). Similar results were 
found for a number of more massive Herbig-Ae/Be stars (e.g., 
|Natta et al.||2004| l. From this mm slope one can estimate the 
opacity index /? (1 -i- A) x (a - 2), where A is found to be 
-0.20 ( |Rodmann et al.||2056l [Lommen et alT||2007l l, and values 
of P ^ 1 .0 for A> \ mm were found. Such a slope can be natu- 
rally explained by a significant fraction of grains at least several 
mm in size present in the discs ( Draine|2006 1. 

A subsample of the sources observed by Lommen et al 



( 2007 [ l overlapped with the Spitzer Infrared Spectrograph (IRS) 



observations pu bUshed by [Kessler-Silacci et "aL] P006 ) and 
Lommen et al.| p007 ) found a tentative correlation between 
the mm slope of the spectral energy distribution (SED) and 
the strength and shape of the 10-yum silicate feature for these 



sources. Note that the 10-/zm feature primarily probes the hot 
surface layers of the inner disc, whereas the (sub)mm obser- 
vations provide information of the cold mid-plane of the outer 
disc. A correlation between the two is therefore not obvious at 
all and a confirmation of this correlation would give very valu- 
able information on the processes of dust growth in protoplane- 
tary discs, as it would imply that grain growth from submicron 
to mm sizes is both fast and occurs simultaneously throughout 
the whole disc. 



Acke et al. (2004i calculated the (sub)mm spectral indices 



of 26 Herbig-Ae/Be stars, for which the infrared SED could also 
be determined. They found a correlation between the strength 
of the ratio of the near- to mid-infrared excess and the slope of 
the (sub)mm energy distribution for these sources, which they at- 
tributed to a correlation between the disc geometry (flared versus 
self-shadowed) and the size of the grains in the disc. However, 
the authors did not find a correlation between the strength and 
the shape of the 10-yum silicate feature and the (sub)mm spectral 
index (see also | Acke & van den Ancker|2004| l. 

The aim of this paper is to investigate the tentative correla- 
tion between the strength and shape of the IQ-jiva silicate feature 
and the spectral slope in the (sub)mm regime, found by Lommen 
et al.| ( |200 7), for a larger sample. A subsample of sources studied 
with the Spitzer IRS were observed with mm and cm interferom- 
eters (Section 2). Interferometers were used to ascertain that the 
emission is dominated by disc emission, since extended emis- 
sion from surrounding material will be filtered out. Also, spa- 
tially resolving the disc ensures that the emission is not optically 
thick (e.g., |Natta et al.||2004|). The results of the observations, 
including dust disc masses and mm slopes, are shown in Section 
3, and in Section 4 we present model results for discs. The ob- 
servations and models are compared and discussed in Section 5; 
conclusions are formulated in Section 6. 



2. Observations 

For this study, we compared Spitzer IRS observations cover- 
ing the 10-/im silicate feature with mm observations from the 
Very Large Array (VLA, operated by NRAcQi, the Combined 
Array for Research in Millimeter-wave Astronomy (CARMA|^, 
the Submillimeter Array (SMA[^, and the Australia Telescope 
Compact Array (ATCAfl). The sources for which new obser- 
vations are obtained for this work are listed in Table |2] A full 
log of the newly obtained mm and cm observations is listed in 
Appendix [a] A full log of the newly obtained mm and cm results 
is listed in Appendix [B] 



' The National Radio Astronomy Observatory is a facility of the 
National Science Foundation operated under cooperative agreement by 
Associated Universities, Inc. 

- Support for CARMA construction was derived from the Gordon 
and Betty Moore Foundation, the Kenneth T. and Eileen L. Norris 
Foundation, the Associates of the California Institute of Technology, the 
states of California, Maryland, and Illinois, and the National Science 
Foundation. Ongoing CARMA development and operations are sup- 
ported by the National Science Foundation under a cooperative agree- 
ment, and by the CARMA partner universities. 

^ The Submillimeter Array is a joint project between the Smithsonian 
Astrophysical Observatory and the Academia Sinica Institute of 
Astronomy and Astrophysics and is funded by the Smithsonian 
Institution and the Academia Sinica. 

"* The Australia Telescope Compact Array is part of the Australia 
Telescope which is funded by the Commonwealth of Australia for op- 
eration as a National Facility managed by CSIRO. 
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Table 1. Distances to and ages of star-forming clouds. 



Cloud 


Age" 


O" 




(Myr) 


(pc) 


Lupus 1 and 2 


< 1 


150 ± 20 


Lupus 3 


1-1.5 


200 ± 20 


Lupus 4 


1-1.5 


165 ± 15 


Chal 


3^ (southern subcluster) 


160 ± 15 




5-6 (northern subcluster) 




Corona Australis 


5-13 


-130 


Serpens 


1-15 


259 ± 37 


Gum nebula 


2-6 


400 ± 60 


Taurus-Auriga 


1-10 


140 ± 15 



Ages adopted from |H eiles ( 1998), Comeron et al. j2003[l, |James 
|et al.| pOOell, |Luhniian| ^2007; , ,Comer6n| (|2006|>, Keiiyon et al. 
( |2008| >, TNeuhauser & Forbrich| ( |200"8| >, |01 veira et al.H2009> , md 

references therein. 

Distances adopted from 



Brandt et al. 



19711, Kenyon et al. i 1994 



Straizys et al.|fT 996^,'Whittet et al. n997VBertout et al. ( 



mm 



de Zeeuvy et al.| ( 19991, Com eron, (2006j , .Neuhauser & Forbr ich, 
( |2008| ), and references therein. 



2. 1. Source selection and Spitzer data 

To look for possible environmental effects, sources in a total of 
five star-forming regions were observed, spread over the constel- 
lations Lupus, Chamaeleon, Corona Australis, Serpens, and the 
Gum nebula in Vela at distances of about 150-200, 160, 130, 
260, and 400 pc, respectively. Furthermore, data from the lit- 
erature for the Taurus-Auriga star-forming region at about 140 
pc were included to improve the statistics further, see Table T] 
The sources were pre-selected to have a large spread in the 
strengths and shapes of the lO-yLim features from Spitzer IRS 
data, mainly the "From Molecular Cores to Planet-forming 
Discs" programme (c2d, ,Evans et aL| 2003 Program IDs 139 
and 172-179), the "The evolution of dust mineralogy in south- 
ern star forming clouds" programme (CM. Wright PI, Project 
ID 20611), and "A complete IRS survey of the evolution of 
circumstellar disks within 3 Myr: New clusters of sequential 
star formation in Serpens" (K.M. Pontoppidan PI, Project ID 
30223). The spectra from the c2d project were previously pub- 
lished in'Kes sler-Silacci et al.| ( |2006l l and |01ofsson et al.| ( |2509j ). 
Program P2061 1 includes Spitzer IRS observations from embed- 
ded YSOs, T-Tauri stars, and Herbig/Vela-type stars. The results 
for the T-Tauri stars are presented in this work. 

The data from Project ID 20611 are presented here for the 
first time. The data from the other programmes are re-reduced 
for this work using the updated c2d IRS reduction pipeline 
( Lahuis et al.|2006) l for uniformity of the comparisons. Spectra 
were obtained both integrated over the full aperture of the in- 
strument as well as convolved with the point spread function 
(PSF) at each wavelength. The spectra obtained using the Full- 
Aperture extraction method were used in here, unless the final 
spectrum quality of the PSF extraction method was consider- 
ably better. Furthermore, only data from the short-low module 
(SL, 5.2-14.5 yum) were included, unless data from the short- 
high module (SH, 9.9-19.6 fim) were present and of significantly 
higher quality. 

In binary systems, it is possible that circumstellar discs get 
truncated due to binary interaction, affecting grain growth in the 
discs. To check for such effects, a number of binaries were in- 
cluded in the sample. Furthermore, the sources were selected 
to include so-called "cold" or "transitional" discs (e.g. |Brown| 



|etal.|2007l l. The cold discs show a lower flux in the mid-infrared, 
which can be naturally explained by a lack of small warm dust 
close to the star. Several of the cold discs were recently found 
to be circumbinary discs, with a large hole or gap in the centre, 
e.g., CS Cha ( [Espaillat et al.|2007| l and HH 30 ( |GuiUoteau etaT 



2008 ). However, some cold discs are supposedly single stars, re- 
quiring a different mechanism to clear the inner discs of small, 
hot grains (e.g., Pontoppidan et al. 2008) . One such mechanism 
could be grain growth into larger particles. Another possibility 
would be that a planet has cleared the inner disc from most of 
the large grains, leaving behind a protoplanetary disc dominated 
by small, micron-sized grains. A number of cold discs of Brown] 
|et al.,(2007) and Merin et al. (2008 1 were included in the sample 
with the aim to explore this possibility. A full list of the sources 
(35 single sources and five binaries) is given in Table 2] As will 
be shown in the next Section, 33 of these turn out to have a de- 
tected 10-jum feature and 13 yield a mm slope, more than dou- 
bling the sample of sources studied in Lommen et al. ( 2007| l. 



2.2. SMA observations 

15 single sources and one binary were observed with the SMA 
for the project 2007B-S033. The observations were carried out 
on 14 Mai-ch and 19 April 2008. The data of 14 March were 
unusable due to phase instabilities and the track was reobserved 
on 7 May 2009. On 19 April 2008, the phases were stable and 
the zenith optical depth at 225 GHz was around T225 =0.13 all 
through the night. The synthesised beam was about 4.8 x 2.8 arc- 
sec (natural weighting). On 7 May 2009, the phases were stable 
and T225 was low with values ranging from 0.05 to 0.08. The 
synthesised beam was about 4.1 x 2.2 arcsec. The two side- 
bands were combined into one continuum channel to improve 
the signal-to-noise ratio, resulting in an effective wavelength of 
1.33 mm. 

The sources VV CrA (binai-y), S CrA (binary), and DG CrA 
(single source) were observed as part of the SMA "filler" project 
2008A-S111 on 1 October 2008. Only six of the eight antennas 
were available for this track. However, T225 ~ 0. 1 and the phases 
were stable, resulting in extremely good data. The synthesised 
beam of the resulting maps was about 5.0 x 2.1 arcsec (natu- 
ral weighting). The correlator was tuned to 218 and 228 GHz. 
Combination of the two sidebands resulted in an effective wave- 
length of 1.35 mm. 

The absolute flux calibration of the first track (19 April 2008) 
was carried out on Mars and the resulting fluxes are estimated to 
be accurate to about 20%. The second and third tracks (1 October 
2008 and 7 May 2009) were flux calibrated on CaUisto. The un- 
certainty in the absolute fluxes for those tracks is estimated to be 
15% or better 

Hence, a total of 16 single sources and three binaries located 
in the Lupus star-forming region were observed with the SMA 
for this project. The sources are listed in Table |2] a detailed log 
of the observations is given in Table A. 1 and detailed results are 
presented in Table B.l and Fig. B.l 



2.3. ATCA observations 

The data for the ATCA project CI 794 were taken over the pe- 
riod July to August 2008 when the array was in the H214 con- 
figuration. A total of 15 sources were observed: 14 sources (in- 
cluding the binary IK LupH-Sz 66) were measured at 3 mm and 
1 1 sources at 7 mm. The sources are listed in Table |2] a de- 
tailed log of the observations is given in Table] A. 2] and detailed 
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Table 2. List of sources observed with the SMA, ATCA, CARMA, and VLA. 



Source 


RA 
(J2000) 


Dec 
(J2000) 


Sp. T. 


Cloud 


Spitzer 


1 mm 


3 mm 


7 mm 


Comments 


Vela 


HBC 553 


08 08 22.2 


-36 03 47.0 


M1.5 


Vela 


Wright 










HBC 556 


08 10 30.9 


-36 01 46.5 


M4 


Vela 


Wright 




ATCA 


ATCA 




HBC 557 


08 12 47.0 


-36 19 18.0 


K3: 


Vela 


Wright 




ATCA 


ATCA 




HBC 559 


08 13 56.1 


-36 08 02.1 




Vela 


Wright 




ATCA 


ATCA 




HBC 560 


08 14 21.9 


-36 10 03.4 


K8 


Vela 


Wright 










HBC 561 


08 15 55.3 


-35 57 58.1 


K8 


Vela 


Wright 








Binary" 


Chamaeleon 


SZ Cha 


10 58 16.9 


-77 17 17.6 


KOe 


Cha I 


GTO 




ATCA 






Sz 32 


11 09 53.4 


-76 34 25.5 


K4.7 


Cha I 


c2d 




ATCA 


ATCA 




Lupus 


IKLup 


15 39 27.8 


-34 46 17.2 


K7 


Lupus 1 


Wright 


SMA 


ATCA 




Binary'' 


Sz 66 


15 39 28.3 


-34 46 18.0 


M2 


Lupus 1 


Wright 


SMA 


ATCA 




Binary'' 


HMLup 


15 47 50.6 


-35 28 35.3 


M4 


Lupus 1 


c2d 


SMA 








Sz 73 


15 47 56.9 


-35 14 34.7 


MO 


Lupus 1 


c2d 


SMA 








HNLup 


15 48 05.2 


-35 15 52.8 


M1.5 


Lupus 1 


Wright 


SMA 








Sz 76 


15 49 30.7 


-35 49 51.4 


Ml 


Lupus 1 


c2d 


SMA 








Sz 77 


15 51 47.0 


-35 56 42.8 


MO 


Lupus 1 


Wright 


SMA 








IMLup 


15 56 09.2 


-37 56 05.9 


MO 


Lupus 2 


c2d 


SMA"^ 


ATCA'^ 


ATCA 




RYLup 


15 59 28.4 


-40 21 51.2 


GOV: 


Lupus 3 


c2d 


SMA 


ATCA 


ATCA 




MYLup 


16 00 44.6 


-41 55 29.6 




Lupus 4 


Wright 


SMA 


ATCA 


ATCA 




EX Lup 


16 03 05.5 


-40 18 25.3 


MO 


Lupus 3 


c2d 


SMA 








Sz91 


16 07 11.6 


-39 03 47.1 


M0.5 


Lupus 3 




SMA 






Cold disc 


Sz 96 


16 08 12.6 


-39 08 33.3 


M1.5 


Lupus 3 


c2d 


SMA 








Sz 102 


16 08 29.7 


-39 03 11.0 


K:e 


Lupus 3 


c2d 


SMA 








Sz 111 


16 08 54.7 


-39 37 43.1 


M1.5 


Lupus 3 




SMA 


ATCA 


ATCA 


Cold disc 


SSTc2d J161029.57-392214.7 


16 10 29.6 


-39 22 14.7 




Lupus 3 


c2d 


SMA 






Cold disc 


SSTc2dJ161 159.81-382338.5 


16 11 59.8 


-38 23 38.5 




Lupus 3 


c2d 


SMA 








RXJ1615. 3-3255 


16 15 20.2 


-32 55 05.0 


K5 


Isolated 


c2d 


SMA 


ATCA 


ATCA 




Serpens 


VV Ser 


18 28 47.9 


-1-00 08 40.0 


A2e 


Serpens 


c2d 


CARMA 


CARMA 


VLA 


Herbig Ae 


SSTc2d J182850.20-^000949.7 


18 28 50.2 


-1-00 09 49.7 


M5+4 


Serpens 


c2d 


CARMA 


CARMA 


VLA 




SSTc2d J182858.08-^001724.4 


18 28 58.1 


-1-00 17 24.4 


G3±5 


Serpens 


c2d 


CARMA 


CARMA 




Cold disc 


SSTc2d J182900.88-H002931.5 


18 29 00.9 


-1-00 29 31.5 


K7±2 


Serpens 


c2d 


CARMA 


CARMA 


VLA 




CoKu Ser-G3 


18 29 01.8 


-1-00 29 54.6 


K0±7 


Serpens 


c2d 


CARMA 


CARMA 


VLA 




IRAS 18268-0025 


18 29 28.1 


-00 22 58 




Serpens 


c2d 


CARMA 


CARMA 


VLA 




SSTc2d Jl 82936. 1 9-^0042 16.7 


18 29 36.2 


-1-00 42 16.7 


F9±5 


Serpens 


c2d 


CARMA 


CARMA 






SSTc2d J182944. 10-^003356. 1 


18 29 44.1 


-1-00 33 56.1 


M0±1.5 


Serpens 


Pontoppidan 


CARMA 


CARMA 




Cold disc 


EC 82 


18 29 56.8 


-1-01 14 46.0 


MO 


Serpens 


c2d 


CARMA 


CARMA 


VLA 




EC 90 


18 29 57.5 


-1-01 14 07 


M4±2 


Serpens 


c2d 


CARMA 


CARMA 


VLA 


Binary'' 


EC 97 


18 29 58.2 


-1-01 15 22 




Serpens 


c2d 


CARMA 




VLA 




GSC 00446-00153 


18 30 06.9 


-1-00 42 34 


F3V 


Serpens 




CARMA 


CARMA 






Corona Australis 


SCrA 


19 01 08.6 


-36 57 20.0 


K3 


CrA 


Wright 


SMA 


ATCA 


ATCA 


Binary"^ 


DG CrA 


19 01 55.2 


-37 23 40.5 




CrA 


Wright 


SMA 


ATCA 


ATCA 




VV CrA 


19 03 06.7 


-37 12 49.7 


K7 


CrA 


Wright 


SMA 


ATCA 


ATCA 


Binary*^ 



Separation a.' 63 ( |Correia et al.|2006| l. 

IK Lup (Sz 65) and Sz 66 form a binary with a separation of 6'.'4 (this work). 
Published in |Lomm en et al.'(20 07|l. 

Separation r/5 (e.g.jciardi et al.^005 and references therein). 
Separation V.'3 (e.g., Forbr ich et al.|2007l and references therein). 
Separation 2'.'0 (this work). 



results are presented in Table B.2 and Figs. B.3 and B.4 The 
weather changed considerably over the course of the observa- 
tions. A short indication of the circumstances for each day is 
included in Appendix [A] Physical baselines ranged from 82 to 
247 metres, resulting in synthesised beam sizes of about 2 arc- 
sec at 3 mm and about 4 arcsec at 7 mm. Combining the two 
sidebands in the 3 mm band resulted in an effective wavelength 
of 3.17 mm, those taken in the 7 mm band in an effective wave- 
length of 6.82 mm. 



The absolute flux calibration for the first track was canied 
out on Mars, whereas the flux calibration for the other tracks was 
carried out on Uranus. Only the shortest baselines were taken 
into account when determining the absolute gain offset so as 
to minimise the possible effect of the planets' being resolved. 
Furthermore, the planets were observed at elevations close to 
those at which the gain calibrators were observed. Overall, the 
uncertainty in the absolute fluxes is estimated to range from 15 
to about 25%. 
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2.4. CARMA observations 

For this work, eleven single sources and one binary located in 
Serpens were observed with CARMA at 1 and 3 mm in the 
period April to June 2008 for project c0165. The sources are 
listed in Table |2j a detailed log of the observations is given in 
Table A. 3 and detailed results are presented in Table B.3 and 
Figs. 



B.5 and IB. 61 Weather conditions varied over the course of 



the observations, with a typical water path length of 3-6 mm. 

The gain calibrator originally selected for the observations 
at 1 mm, QSO J1743-038, turned out to be too weak to perform 
a decent gain calibration, rendering most of the C-configuration 
observations unusable. For the second part of the observations 
the telescope was in the D configuration (baselines 11-148 me- 
tres), yielding a synthesised beam of about 3x2 arcsec at 1 mm 
and about 6x4 arcsec at 3 mm. The effective wavelength of the 
1 mm-band observations was 1.33 mm, that of the 3 mm-band 
observations 3.15 mm. 

The absolute fluxes were calibrated on the quasars QSO 
J2253H-161 (3c454.3), QSO J1229H-020 (3c273), and QSO 
J 1256-057 (3c279), whose fluxes were bootstrapped from planet 
observations on short baselines on dates as close as possible to 
the observation dates. The fluxes of these quasars vary consid- 
erably over the course of weeks to months at 1 and 3 mm, but 
day-to-day variations are usually less than 10%. Taking this into 
account, the effective uncertainty in the absolute fluxes for our 
target sources is estimated to be less than 30%. 

2.5. VLA observations 

Of the sources in the Serpens star-forming region observed with 
CARMA, seven single sources and the binary EC 90 were ob- 
served with the VLA at 7 mm and at 1.3, 3.6, and 6.3 cm under 
programme AL720. The sources are listed in Table [2] a detailed 
log of the observations is given in Table A.4 and detailed results 



are presented in Table B.4 The observations were carried out 



from 10-15 March 2008, when the array was in the C configu- 
ration, with baselines of up to 3.6 km and a synthesised beam of 
about 0.5 arcsec at 7 mm. All observations were performed in 
the default continuum mode in which, at each frequency, the full 
100-MHz bandwidth was used in two adjacent 50 MHz bands. 
Although weather conditions were good in general, a few hours 
of observing time were lost at the end of the last two tracks due 
to high winds. 

The VLA data were flux calibrated on the quasar QSO 
J1331H-305 (3c286). The flux as a function of wavelength is 
modelled by the AIPS reduction package. The resulting uncer- 
tainty in the absolute flux calibration is estimated to be about 
20% at 7 mm and 1.3 cm and better than 10% at 3.6 and 6.3 cm. 



3. Results 

3. 1. Mm and cm source fluxes and dust disc masses 

A full log of the results is listed in Appendix B. The results of 
the interferometric observations at 1, 3, and 7 mm are listed in 
Table [3] 

A total of 16 single sources and three binaries in Lupus are 
observed with the SMA. Nine of the single sources are detected 
and one of those, Sz 73, turned out to harbour two sources, with 
a projected separation of about 4 arcsec. It is possible that the 
detection of Sz 73 with SEST (Nurnberger et al.|1997| l included 
both sources. The binaries VV CrA and S CrA are detected and 
unresolved. Of the binary system IK Lup (Sz 65) and Sz 66, only 



IK Lup is detected, although a second peak is detected at 2 arc- 
sec from the 2MASS position of Sz 66. Sz 66 was previously 
detected with a S /N of almost four using the SEST bolometer. 
All sources in Lupus observed with the ATCA at 3.2 mm are de- 
tected; the binary system IK Lup and Sz 66 remained unresolved. 
Only one Lupus source, IM Lup, is detected at 6.8 mm. MY Lup 
would have been detected at 6.8 mm with a signal-to-noise ratio 
of about ten if it had a similar mm slope as IM Lup. 

None of the three sources in the Gum nebula observed with 
the ATCA at 3.2 and at 6.8 mm are detected at either wavelength 
down to 3cr upper limits of ~3 mJy at 3.2 mm and of ~0.5 mJy 
at 6.8 mm. This can be attributed to the large distance between 
us and this star-forming region. If the sources in the Gum neb- 
ula had similar luminosities as those in the Lupus clouds, they 
would have had a flux of ~0.7 mJy at 3.2 mm, which is below 
the noise level. Note that, although the Vela molecular ridge has 
been observed at mm wavelengths ( jMassi et al.|1999| |2()()7[ |. no 
published mm continuum data of the Gum nebula exist in the 
literature. 

The source SZ Cha is detected at 2.3 mJy at 3.2 mm. Sz 32 
is not detected down to a 3cr upper limit of 2.9 mJy at 3.2 mm. 
It is, however, detected with a flux of 0.77 mJy at 6.8 mm. 

VV CrA and S CrA are clearly detected at 1.3 mm with the 
SMA, with fluxes of 376 and 303 mJy. DG CrA, however, is not 
detected, down to a 3cr upper limit of only 6.6 mJy. VV CrA and 
S CrA are also easily detected with the ATCA at 3 and 7 mm. 

Of the sources in the Serpens star-forming region that were 
observed with the CARMA, only three are detected: the single 
sources SSTc2d J182900.88H-002931.5 and GSC 00446-00153 
and the binary system EC 90, which remained unresolved. This 
can in part be explained by the distance to the star-forming 
region in Serpens, which is larger than those in Chamaeleon, 
Lupus, and Corona Australis. Furthermore, some of the sources, 
of which six are new Spitzer sources, may have an intrinsically 
lower luminosity. None of the sources are detected at 6.8 mm 
using the VLA. 

Four cold discs are observed at 1.3 and 3.2 mm for this work. 
Only one of those, Sz 1 1 1, is detected. Unfortunately, Sz 1 1 1 was 
not observed with the Spitzer IRS. 

All four binaries that are observed at 1.3 and 3.2 mm are 
detected at both wavelengths. However, in the case of the bi- 
nary consisting of the stars IK Lup and Sz 66, only the for- 
mer is detected at 1.3 mm. EC 90 and S CrA remain unre- 
solved. The binary IK LupH-Sz 66 is resolved with the ATCA 
at 3.2 mm. VV CrA is resolved with a binary separation of 
2.0 arcsec with the ATCA at 3.2 mm if the source is imaged 
using uniform weighting (optimised for resolution). However, 
this binary remains unresolved with the SMA at 1.3 mm (beam 
size 4.7 X 1.9 arcsec) and with the ATCA at 6.8 mm (beam size 
5.3 X 3.0 arcsec). 

The detection rate of the sources observed in this study is 
rather low. This can in part be understood by the distance to the 
star-forming regions, with the Serpens star-forming region being 
almost twice as far away as the Taurus-Auriga star-forming re- 
gion and the Gum nebula in Vela almost three times as far away. 
This reduces the observed flux for similar sources by a factor of 
about four to nine. The low detection rate for Lupus is largely a 
selection effect: most of the brightest sources had been observed 
before ( Lommen et al.|2007 1. These previously detected sources 
were not reobserved for this work, but their published values will 
be included in the analysis below. 

Dust disc masses are obtained from the fluxes at 3.2 mm un- 
der the rather crude assumptions of an isothermal disc and a fixed 
opacity. Assuming also optically thin mm emission, the dust disc 
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Table 3. Fluxes from point-source fits in the (m, v) plane obtained from interferometric data and single-dish 1.20-1.27 mm SEST 
fluxes. 



Source 


1.3 mm" 


3.2 mm" 


D.C 


mm 


SEST 1.20-1.27 mm*-' 




Flux 


rms 


Flux 


rms 


Flux 


rms 


Flux 


rms 




(mJy) 


(mJy/bm) 


(mJy) 


(mJy/bm) 


(mjy; 


(mJy/bm) 


(mJy) 


(mJy) 


HBC 556 


— 


— 


< 3.7" 


1.2 


< 0.7" 


0.22 


— 


— 


HBC 557 


— 


— 


< 3.2" 


1.1 


< 0.6" 


0.18 


— 


— 


HBC 559 


— 


— 


< 2.9" 


1.0 


< 0.3" 


0.09 


— 


— 


SZ Cha 


— 


— 


2.3 


0.4 






77.5 


20.3 


Sz 32 


— 


— 


< 2.9" 


1.0 


0.77 


0.14 


93.1 


20.8 


IK Lup 


28 


2.8 


3.4 


0.4 






56 


10 


Sz 66 


< 8" 


2.8 


2.2 


0.4 






47 


12 


HM Lup 


< 10" 


3.4 


— 


— 






<45" 


15 


Sz 73 a 


16 


2.9 


— 


— 






26 


8 


Sz 73 b 


16 


2.9 


— 


— 






— 


— 


HN Lup 


15 


3.0 


— 


— 






< 51" 


17 


Sz 76 


12 


3.3 


— 


— 






< 45" 


15 


Sz 77 


< 10" 


3.2 


— 


— 






< 45" 


15 


iM Lup 


188 


4.3 


8.9 


1.3 


2.2 


0.16 


260 


9 


RY Lup 


78 


4.9 


2.8 


0.7 


< 0.6" 


0.21 






AyTV T nr» 
iVl I LUp 


56 


3.4 


8.7 


0.4 


< u.o 


u.zu 


— 


— 


mp 


19 


3.9 














Q-7 01 


< 13" 


4.3 










< 27" 


9 


oZ yO 


< 13" 


4.2 










< 45" 


15 


oZ iUZ 


< 11" 


3.8 










< 30" 


10 


^■7 1 1 1 

oZ 111 


49 


4.8 


5.7 


0.7 


< u.o 


U. 1 V 








< 13" 


4.4 














oo 1 CZa J iOl Vjy .o\.-JoZ.jjo.O 


< 13" 


4.2 














P V T1 Al ^ 


132 


3.9 


6.7 


0.6 


< U.J 


n 1 7 

U. 1 / 






LC oZ 


< 15.7" 


5.2 


< 2.9" 


LO 


< U.J 


r\ 7 
U./ 


— 


— 


Lu yu 


95.7 


10.4 


11.5 


L2 


< l.U 


U.J 


— 


— 


LC. V / 


< 23.3" 


7.8 






< U.O 


r\ 7 
U./ 






oo ICZtl J ioZ7UU.OO+UUZ!7jl.Z) 


26.4 


4.9 


3.4 


0.5 


< U. / 


n 7 


— 


— 




< 15.6" 


5.2 


< 1.9" 


0.6 






— 


— 


CoKu Ser-G3 


< 17.1" 


5.7 


< 1.8" 


0.6 


< 1.2 


0.4 






SSTc2d J182858. 08+001724.4 


< 24.6" 


8.2 


< 1.9" 


0.6 










VV Ser 


< 14.8" 


4.9 


< 1.8" 


0.6 


< 0.7" 


0.2 






SSTc2d J 182850.20+000949.7 


< 22.8" 


7.6 


< 1.9" 


0.6 


<0.6" 


0.2 






SSTc2d J 1 82944. 10+003356. 1 


< 15.0" 


5.0 


< 1.7" 


0.6 










SSTc2d J 1 82936. 19+0042 16.7 


< 8.7" 


2.9 


< 2.9" 


1.0 










GSC 00446-00153 


90.8 


3.5 


6.8 


1.0 










VV CrA 


376 


4.5 


26.8 


1.1 


8.3 


0.25 


469, 584 


21 


S CrA 


303 


3.2 


24.9 


1.0 


3.7 


0.20 


290 




DO CrA 


< 6.6 


2.2 


< 2.5" 


0.8 


<0.5" 


0.16 







1 mm band observations are at 1.33 mm (SMA, Lupus and CARMA, Serpens) and at 1.35 mm (SMA, Corona Australis). 
3 mm band observations are at 3.17 mm (ATCA, Lupus and Corona Australis) and at 3.15 mm (CARMA, Serpens). 

SEST fluxes are from |Numberger et al.| ( |1997 Lupus), [Henning et al.| ( |199?i Cham aeleon), and,Henning et aL] ( |1994[ Corona Australis), with 
an adopted centre frequency of 236 GHz (1.27 mm). The values in italic are from |Chini et al.| ( |2003^ , with an adopted centre frequency of 
250 GHz (1.20 mm). 
Quoted value is 3(T upper limit. 



mass is given by Mjisc - FyD^ I KyByiJi^si)^ where D is the dis- 
tance to the source, Kv the dust opacity (taken to be 0.9 cm^ g"', 
cf. |Beckwith et al.|1990 ), and Bv(7^dust) the brightness at the fre- 
quency V = 94 GHz for a dust temperature r^ust^ as given by the 
Planck function. We assume a dust temperature Tdust = 25 K and 
find dust disc masses ranging from ~0.4 to ~ 7 x 10""* Mq. The 
dust disc masses are presented in Table |4] 

3.2. Millimetre slopes 

The fluxes at 1, 3, and 7 mm can be combined to obtain the spec- 
tral index a, where Fy oc v". We are interested in the emission 
coming from the dusty disc. However, at 7 mm, other emission 



mechanisms may contribute significantly to the flux. Sources 
may include an ionised wind or chromospheric magnetic activ- 
ity. |RodmannetaLj(2006 ) compare their fluxes at 7 mm to those 
at 3 and 6 cm and claim that about 20% of the emission at 7 mm 
is due to free-free emission. On the other hand, [Lommen et aL] 
( |2009| l find that the emission at 7 mm can be entirely attributed 
to dust emission for a small sample of three sources. It is possi- 
ble that the emission due to, e.g., an ionised wind, is independent 
of the disc mass and thus the relative contribution from such a 
wind will be larger for young stellar objects that are weaker at 
mm wavelengths. This could explain the findings of |Lomrnen 
iet al.| p009) l, who monitored some of the strongest pre-main- 
sequence mm emitters in the southern sky. However, a larger 
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Source 


a 1-3 


Q'3_7 


Dust disc mass" 




peak 










(lO-^Mo) 






HBC 553 








0.22 


1.43 


0.86 


HBC 556 






< 5.1 


0.40 


1.64 


0.99 


HBC 557 






< 4.4 


0.21 


1.37 


0.95 


HBC 559 






< 4.0 


0.39 


1.60 


0.85 


HBC 560 






... 


0.29 


1.44 


1.03 


HBC 561 








0.43 


1.67 


0.90 


SZ Cha 


3.8 ±0.4" 




0.6 


0.71 


2.30 


0.85 


Sz 32 


> 3.7'' 


< 1.8 


< 0.8 


0.14 


1.24 


1.11 


IK Lup 


2.7 ±0.3" 




0.7 


0.36 


1.51 


0.93 


Sz 66 


1.9 ±0.2'' 




0.4 


0.40 


1.62 


0.92 


HMLup 








0.50 


1.81 


0.85 


Sz 73 








0.26 


1.38 


0.95 


HNLup 






... 


0.18 


1.29 


1.03 


Sz 76 






... 


0.45 


1.77 


0.77 


Sz 77 








0.43 


1.63 


0.94 


IMLup 


3.6 ±0.4" 


1.8 ±0.3 


1.7 


0.52 


1.81 


0.85 


RYLup 


3.8 ±0.5 


> 2.0 


0.5 


1.10 


3.16 


0.66 


MY Lup 


2.1 ±0.5 


> 3.5 


1.7 


0.54 


1.87 


0.79 


EX Lup 








0.54 


2.01 


0.73 


Sz 96 








0.78 


2.30 


0.82 


Sz 102 








0.80 


2.36 


0.82 


Sz 111 


2.5 ±0.4 


> 2.9 


1.1 




... 




SSTc2d J161029.57-392214.7'' 






<0.3 


0.29 


1.52 


0.89 


SSTc2dJ161 159.81-382338.5 






<0.3 


0.92 


2.48 


0.84 


RX J1615.3-3255 


3.4 ±0.5 


> 3.4 


1.3 


1.25 


3.04 


0.96 


VV Ser 






< 1.0 


0.33 


1.45 


0.97 


SSTc2d J182850.20+000949.7 






< 1.1 


0.36 


1.57 


0.88 


SSTc2d J182900.88+002931.5 


2.4 ± 0.4 




2.0 


0.15 


1.27 


0.96 


CoKu Ser-G3 






< 1.0 


-0.09 


0.82 


1.17 


IRAS 18268-0025 






<0.7 


0.48 


1.67 


1.14 


SSTc2d 1182944.10-1-003356.1' 






< 0.7 


0.30 


1.53 


0.96 


EC 82 aka CK 3 






< 1.7 


1.41 


3.55 


0.67 


EC 90 aka CK 1 


2.5 ±0.5 




6.7 


0.26 


1.40 


1.12 


EC 97 aka CK 4 








0.34 


1.56 


0.90 


GSC 00446-00153 


3.0 ±0.5 




4.0 








S CrA 


2.9 ± 0.7" 


2.5 ±0.4 


3.6 


0.34 


1.51 


0.97 


DG CrA 






< 0.4 


1.04 


2.94 


0.67 


VV CrA 


2.5 ± oy-^-' ii 


2.4 ± 0.5s 


3.9 









0.9 cm- g"' (cf.|Beckwith et al.|l990|, and a dust 



Dust disc masses estimated from ATCA and CARMA fluxes, assuming a dust opacity k,, ■ 
temperature T^dust ~ 

25 K. 

Using the SEST 1.27 mm (single-dish) flux from'Niimberger et al.Vl 99^ . 
Using the SEST 1.27 mm (single-dish) flux from Henning et al. ( 1993Jl. 
Using the SEST 1.20 mm (single-dish) flux from Chini et al. (2003). 
Using the SEST 1.27 mm (single-dish) flux from Henning et al. ( 1994) ). 
Taking a combined flux of 50.2 mjy at 3. 17 mm. 
Using the PSF extraction method for the Spitzer IRS spectrum. 

Observed with Spitzer IRS for Program ID 30843 (B. Merm PI). The full IRS spectrum will be presented in Merm et al. (2010, in prep.). 
Observed with Spitzer IRS for Program ID 30223 (K.M. Pontoppidan PI). The full IRS spectrum will be presented in Oliveira et al. (2010, 
prep.). 



and more sensitive survey at mm to cm wavelengths is required 
before more quantitative statements on this subject can be made. 
Since we do not have fluxes at all three wavelengths for most 
sources, separate indices will be obtained between 1 and 3 mm 
and between 3 and 7 mm. The results are given in Table |4] 

The slopes between 1 and 3 mm lie between 2.38 ± 0.36 
and 3.83 + 0.46. The opacity index /? can be calculated from 
the mm slope a through /3 ~ (1 + A) x (o- - 2), where A is 
the ratio of optically thick to optically thin emission ([Beckwith 
eraL|[T990j IRodmann et"aL][2506l l. [Rodmann et ai:| ( |2006| l and 
Lommen et al. ( 2007 1 found values of A 0.2 for the sources 



in their samples. Adopting this value, opacity indices /3 of about 
0.5 to 2.2 are found here. The Kolmogorov-Smirnov test gives 
a probability of 50% that the values from this sample and that 
of |Lommen et al.| ( |2007| l are drawn from the same distribution. 
This rather low value is due to the steep slopes for the sources 
RY Lup (3.83 ± 0.46) and SZ Cha (3.78 ± 0.43). Note that the 
corresponding values for /3 are > 2, whereas the value for the 
interstellar medium is ySisM ~ 2.0 (e.g., Draine & Lee|1984 l. 



A mm slope between 3 and 7 mm could only be deter- 
mined for three sources, whereas lower limits are found for 
four more sources and a strict upper limit of ai-^ < 1.8 for 
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Sz 32. Interestingly, a lower limit of ai-j > 3.7 is found for 
Sz 32 between 1 and 3 mm. Other emission mechanisms (due to, 
e.g., a wind or chromospheric activity) may contribute at 7 mm. 
Although it is found that for most sources this contribution is 
only of the order of 20% ( jRodmann et al.| 20061), it is possible 
that it is higher for Sz 32, causing the very shallow slope between 
3 and 7 mm. The slopes of a^-j = 2.4 + 0.5 and 2.5 + 0.4 for 
VV CrA and S CrA are consistent with those of o-i-s = 2.5 ± 0.5 
and 2.9 + 0.7 and also the slopes between 3 and 7 mm found 
for RY Lup (> 2.0), Sz 111 (> 2.9), RX J1615.3-3255 (> 3.4), 
and MY Lup (> 3.5) are consistent with the values between 1 
and 3 mm. The slope between 3 and 7 mm for IM Lup, how- 
ever, is very shallow compared to that between 1 and 3 mm: 
as-T = 1.8 + 0.3 vs ai^3 = 3.6 + 0.4. |Pinte et aL] ( |2008l l found a 
mm spectral index of 2.80 ±0.25 and their modelling results sug- 
gested that IM Lup has grains of at least mm sizes in the disc. A 
shallowing of the slope beyond 3 mm may indicate the presence 
of at least cm-sized grains. A similar effect on the cm SED was 
found for TW Hya ( |Wilner et al.|2005] l. 



3.3. Results from Spitzer infrared observations 

The spectra of the T-Tauri stars observed for Spitzer project 
P20611, including sources in Lupus, Corona Australis, and 
the Gum nebula, are published for the first time here and 
shown in Fig. [T] The spectrum of VV CrA is saturated be- 
low 10 /urn and excluded from the sample. The spectrum of 
SSTc2d J161029.57-392214.7 (P30843, B. Men'n PI) will be 
published in Merm et al., (2010, in prep.). The spectrum of 
SSTc2d J182944.11H-003356.1 (P30223, K.M. Pontoppidan PI) 
will be published in Oliveka et al. (2010, in prep.). 

The 10-jum silicate features were analysed in the ways of 
both'Furlan e t al.| ( |2006l l and |Kessler-Silacci et al.| ( |2006l l. |FurIanl 
[et al. (2006) fitted a third-order polynomial to the continuum 
around the 10-yum feature and determined the integrated flux 
above and below the continuum. The strength of the IQ-fim fea- 
ture was then defined as the ratio of the integrated flux due to 
the feature divided by the integrated flux due to the continuum, 
(Fio - /^cont)/^^cont, resulting in a strength larger than for a fea- 
ture in emission. Kessler- Silacci et al.| ([2006 ) determined the 
continuum in three different ways, depending on the full mid- 
infrared SED and which data were available for each source, and 
subsequently determined the normalised spectra S y according to 



5, = 1 -H 



(1) 



where Fy is the observed spectrum, Fy^c is the fitted contin- 
uum, and < Fyc > is the frequency-averaged continuum flux 
(see their paper for details). They defined the strength of the 10- 
fim feature as the maximum value of Sy between 5 and 13 jjm. 



resulting in a value larger than 1 for a feature in emission. 



I lOfim 

peak 

Furthermore, [Kessler-Silacci et al.[ ( |2006| l defined the shape of 
the 10-yum feature as the ratio 5 11.3/5 g.g- 

For this work, the continuum was consistently chosen for 
all sources by fitting a third-order polynomial to data between 
5.0 and 7.5 fim and between 13.0 and 16.0 fim (cf. Furlan 
[et al.|[2006] l. The regular continuum was used rather than the 
frequency-averaged continuum, resulting in the peak strength 
^]ieak shape Fii 3/F9 g. This does not change the results 

significantly (see [ Kessler-Silacci et al.|2()06 l 



The results are listed in Table [4[ and shown in Fig. [2| The 
upper panel of Fig.jijgives (Fio - Fcont)/fcont vs showing 
a clear correlation between the two definitions of the strength of 




10 //m 



peak 



Fig. 2. Upper panel: the peak strength of the 10-yum feature as 
de fined in[Kessler-Silacci et al.[ ( [2006| l vs the strength as defined 



m 



Furlan et al. 



^(2006). Lower panel: the peak strength of the 

10-;um feature vs the shape as in Kessler-Silacci et aL] ( 2006) . 
Included are the sources observed in this work for which the 
10-/im feature could be obtained. The object towards the top- 
right off the general trend in the lower panel is the Lupus source 
RX J1615. 3-3255, which is isolated from the Lupus clouds. 



the 10-//m feature. The lower panel gives F^^^^ vs Fn^^/Fg s, 
also showing a correlation, confirming the results of KesslerJ 
Silacci et al. ( 2006) ). This figure also confirms that our sample 
covers a large range in silicate-feature characteristics. It follows 
that the three methods to quantify the strength or shape of the 
10-;um feature give completely consistent results. When com- 
paring the 10-/im feature with the mm slope in Sect.[5]below, the 
strength defined as (Fio - f'cont)/^cont will be used. The source 
that lies towards the top and to the right of the general trend in 
the lower panel is RX J1615. 3-3255, an isolated source slightly 
to the north of the Lupus star-forming clouds. 

3.4. lO-jjm feature vs mm slope 

Figs. [3] and [4| show the mm slope a, measured between 1 
and 3 mm, as a function of the strength of the 10-//m feature 
{Fio - /^cont)/^cont- Only the slope between I and 3 mm is used, 
to make the sample as consistent as possible. However, as noted 
before, the slope between 3 and 7 mm is consistent with that be- 
tween 1 and 3 mm for most sources. Included are the sources 
from this study, as well as eleven sources located in the Taurus- 
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Fig. 1. Spitzer IRS spectra from the T-Tauri stars observed for Spitzer project P2061 1 (CM. Wright PI). Spectra with a maximum 
flux below 0. 1 Jy were binned four times to improve the signal-to-noise ratio. 

Table 5. List of sources used in the analysis. 



Source 


(^10 - ^coiil)/^cont 


ffl-3 


Notes 


This work 


RYLup 


1.10 


3.83 




RXJ1615. 3-3255 


1.25 


3.43 




IKLup 


0.36 


3.01 


Binary with Sz 66 


Sz 66 


0.40 


3.29 


Binary with IK Lup 


IMLup 


0.52 


3.60 




SZ Cha 


0.71 


3.78 




SCrA 


0.34 


2.93 


Binary 


EC 90 


0.26 


2.46 


Binary 


SSTc2d J182900.88-^002931.5 0.15 


2.38 




Values taken from^Lommen et al. ( 2007) 


CR Cha 


0.% 


3720"" 




WW Cha 


0.77 


2.70 




HT Lup 


0.29 


2.33 


Binary 


GWLup 


0.40 


2.42 




RULup 


0.39 


2.67 




TCha 


0.06 


2.90 


Cold disc 


CS Cha 


0.91 


2.90 


Cold disc 


DI Cha 


0.39 


> 2.40 




Glass I 


1.02 


> 3.40 






Sources from|Rodmann et al. 


(2006) 




DGTau 


0.00 


2.14 




DOTau 


0.18 


2.29 





Sources from Andrews & Williams ( 2007 



AATau 


0.35 


3.15 


CITau 


0.54 


2.18 


DLTau 


0.09 


1.97 


DMTau 


0.96 


2.91 


DNTau 


0.22 


2.34 


DRTau 


0.27 


2.20 


FTTau 


0.43 


1.82 


GM Aur 


1.19 


3.16 


RYTau 


1.36 


1.91 


AS 205 


0.27 


0.67 



Auriga star-forming region discussed in 'Rodmann et al. ( 2006| l 
and |Andrews & Williams ( [2007j for which the spectral slope be- 
tween 1 and 3 mm could be determined, and the sources located 
in Lupus and Chamaeleon discussed in Lommen et al. ( 2007| l. 
The total number of sources used is 3 1 ; the complete hst is given 
in Table |5] In Fig.|4] the sources are sorted by their star-forming 
region. The smaller symbols designate single stars and the larger 



symbols binaries (or stars that are members of a multiple sys- 
tem). The open symbol to the left is T Cha, an evolved system 
that does not show any silicate emission and is not used in the 
analysis, and the open symbol in the centre designates the "cold 
disc" CS Cha. 

The sources in the sample shown in Fig.|3]lie in a broad band 
roughly running from the lower left (shallow mm slope and weak 
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Fig. 3. The mm slope as measured between 1 and 3 mm as a 
function of the strength of the 10-yum feature. The open sym- 
bols are for T Cha (to the left), which does not show any silicate 
emission and is not used in the analysis, and for CS Cha, a cir- 
cumbinary disc. The dashed line shows a linear fit to all the data. 
The dotted line shows a linear fit to the data with RY Tau ex- 
cluded. Included are the sources studied in this work, as well as 
112006 ), Andrews & WilHams (2007), 



those from Rodmann et al. 



and Lommen et al. (2007 



typical uncertainties. 
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Fig. 4. Same as Fig. [s] with the different sources sorted by 
star-forming region: filled circles: Lupus, five-pointed stars: 
Chamaeleon, cross: Corona Australis, diamonds: Taurus -Auriga, 
and squares: Serpens. The ellipses show the concentrations of 
sources located in the Taurus-Auriga star-forming region (lower 
left), the Chamaeleon 1 cloud (top centre), and the Lupus 1 and 
Lupus 2 clouds (upper left). The remaining two Lupus sources in 
the upper right are an isolated source (RX J1615-. 3-3255, right- 
most dot) and a source from the Lupus 3 cloud (RY Lup, upper- 
most dot). The small symbols designate the single stars and the 
large symbols designate multiple systems. The open five-pointed 
star to the left is for T Cha, an evolved cold disc which shows 
no silicate emission around 10 jim. The open five-pointed star in 
the centre is CS Cha, a circumbinary disc. 



lO-fim feature) to the upper right (steep mm slope and strong 
10-;um feature). The sole exception is the source RY Tau, which 
lies in the lower right corner. The mm slope and the strength 
of the 10-yum feature correlate weakly for the full sample: the 
Spearman rank correlation coefficient is 0.50, with a 99.5% con- 
fidence level. However, if the point for RY Tau is excluded, the 
Spearman rank coefficient becomes 0.66, with a 99.99% confi- 
dence level. Note that RY Tau is a peculiar source: it is found 
to be a rapidly rotating UX Or-type star powering a microjet 
(e.g., |Petrov et al.| 1999 [| Agra- Amboage et al.|20"09| . A possible 
explanation for its location in the 10-yum-feature vs mm-slope 
diagram is a rather evolved disc in which a recent collision event 
produced small grains. This may be similar to the effect recently 
observed in EX Lup, in which a significantly more crystalhne 
lO-jjm feature was observed after an outburst (Abraham et al. 
|2009 ). RY Tau will not be included in the further discussion. 

Fig. |4] suggests a grouping in the /im-vs-mm diagram ac- 
cording to parental cloud, with the sources from the Taurus- 
Auriga star-forming region more concentrated in the lower left, 
the Lupus sources more to the upper left, and the Chamaeleon 
sources more to the centre right. Note that the six Lupus sources 
that are on the left part of the diagram (from top to bottom: 
IM Lup, Sz 66, Sz 65, RU Lup, GW Lup, and HT Lup) are 
all located in the Lupus 1 and Lupus 2 clouds, whereas the re- 
maining two Lupus sources are located in Lupus 3 (top-most 
source, RY Lup) and off-cloud (RX J1615.3-3255). Larger- 
number statistics are needed to confirm this grouping by star- 
forming region in the /im-vs-mm diagram. 

[Kessler-Silacci et aL] ( |2006| found a correlation between the 
spectral type of a source and the strength and shape of the 10-jum 
silicate feature, brown dwarfs having predominantly flatter and 
Herbig-Ae/Be stars having more peaked features. It was found 
that this is most likely due to the location of the silicate emis- 
sion region: [Kessler-Silacci et alT| (|2007) showed that the radius 
of the lO-yum siUcate emission zone in the disc goes roughly as 
{Lt I Lof -^^ . Hence, the feature probes a radius further 

from the star for early-type stars than for late-type stars. In this 
context it is interesting to see whether a correlation with spectral 
type is found in the 10-;um-feature vs mm-slope diagram (Fig.|5]l. 
The M stars in the sample presented here are concentrated to the 
left, the F and G stars to the lower left, and the K stars are found 
both in the lower left and in the upper right. Hence, no clear 
correlation with spectral type is found here. It is interesting to 
note, though, that the F and G sources RY Tau and RY Lup show 
up isolated from the other F and G sources. This may indicate 
that these sources are indeed different from the other sources in 
the sample, justifying the choice not to include RY Tau in the 
analysis. 



4. Modelling 

4. 1. Disk model parameters and SEDs 

Variations in the strength and shape of the IQ-fim feature (e.g., 
2006| l as well as in the (sub)mm slope (e.g.. 



Kessler-Silacci et al. 



Beckwith et al. 1990 1 can be explained by variations in the dom- 
inating grain size in the circumstellar discs, so that one may ex- 
pect a correlation between properties of the 10-;um feature and 
the mm slope. Such a correlation is found for the sample as 
a whole (see previous Section) and this may imply that grain 
growth occurs in the whole disc simultaneously, or that grains 
grow in the inner disc and the new grain size distribution is very 
efficiently spread to the outer disc through radial mixing. Both 
processes will have the effect of a shift of dust mass from small 
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Fig. 5. Same as Fig. [s] with the different sources sorted by spec- 
tral type: circles: F and G, five-pointed stars: K, and squares: 
M. 



particles to larger grains. To study this more quantitatively, we 
ran a number of models with varying grain size distributions. 

We use the axisymmetric radiative-transfer code RADMC, 
developed by Dullemond (20021 and Dullemond & Dominik 
(12004) 



The model consists of a flaring disc, heated passively 
by radiation from the central star, and includes a hot inner wall, 
which is directly irradiated by the central star (Natta et al. 2001, 
see also Dullemond et al. 2001). The surface density of the disc 
as a function of radius Z(r) is defined to be: 



Y(r)^^,,,(r/R,,,r, 



(2) 



with n - -I. The total gasH-dust disc mass was fixed to Mdisc - 
5 X 10"^ Mq. The gas-to-dust mass ratio ^ was set to 100, which 
implies that the total mass in dust grains is also fixed, to a value 
of 5 X 10"^ Mq. The inner radius is fixed to =0.1 AU and 
the outer radius Rom is varied as outlined below. The photons, 
originating from the central star, are allowed to move in three 
dimensions within the axisymmetric grid. In the models, the ge- 
ometry of the disc is determined by vertical thermal equilibrium. 
The temperature and luminosity of the central source are varied, 
as are the power-law slope of the grain size distribution, the disc 
outer radius, and the maximum grain size. 

For the dust opacities, we use a mixture of 80% amorphous 
olivine and 20% armorphous carbon (percentages by mass). The 
opacities are calculated using a Distribution of Hollow Spheres 
(DHS, see |Min et al.|2003] l. The total volume of the spheres oc- 
cupied by the inclusion / is taken in the range / = [0, 0.8] . It was 
found (e. g.'Chiang et al.|2001[|FAlessio et al.|2001[|Dullemond 
|& D ominik 2004) that the mm slope changes if one goes from 
a disc with only "small" particles to a disc that also contains 
some "large" grains. Dullemond & DominikI started with a disc 



in which the dust is made up of only 0.1-jum-sized particles and 
subsequently replaced 90%, 99%, 99.9%, 99.99%, and 99.999% 
of the dust by large, 2-mm-sized grains. The mm slope changes 
considerably when the mass fraction in large grains is changed 
from to 90%, but it does not change further if a larger fraction 



of the dust mass is put in large grains (see Fig. 7 in DuUemond 
|& Domimk||2004[ ). This is a result of the fact that at 1 mm the 
opacity is dominated by the large grains, virtually independent 
of the mass fraction (K. Pontoppidan, priv. comm.). Although it 
is possible that a more gradual change in the mm slope is seen 



Parameter 


standard 


range 


Mass Msiai- 


1.0 Mo 






3000 K 


3000 K, 4000 K 


Luminosity Ls,a,. 


1.0 Lo 


1.0 Lo, 6.0 Lo 


Radius i?stai 


3.7 Ro 


3.7 Rg, 5.1 Lo 


Mass Mjisc (Mq) 


5 X 10-' 




Disc inclination angle / (°) 


45 


5.7, 15, 30, 45, 60, 75 


Surface mass density gradient n 


-1.0 




Gas-to-dust ratio T 


100 




Inner radius Ri„ (AU) 


0.1 




Outer radius ^out (AU) 


300 


100, 200, 300 


Minimum grain size (jim) 


0.003 




Maximum grain size a,„ax (y"m) 




0.1, 1, 10, 100, 






1000, 10000 


Power-law slope m 


3.5 


2.5, 3.0, 3.5, 4.0 



when smaller mass fractions are put in large grains, it does seem 
to be more important what the largest grain size is, rather than 
which fraction of the dust is contained in such large grains. We 
therefore chose not to use a bimodal dust distribution, but a dis- 
tribution in which the size of the grains ranges from a minimum 
value fl^in to a maximum value Omax according to 



n(a) I 



(3) 



This power-law distribution is expected on theoretical grounds 



whenever grain-grain collisions lead to shattering (Dohnanyi 



1969 1. It should be noted that models which include grain growth 
may lead to different grain size distributions (e.g., [Dullemond &| 
Dominik 2005, Tanaka et al. 2005). The value was fixed to 
0.003 fini. Note that ^/min can still have a considerable effect on 
the IQ-fim feature, with large values for amin giving significantly 
weaker features Olofsson et al. (2009). This effect is strongest 
for steep grain size distributions, e.g., m - 4.0. The maximum 
grain size Omax was varied in steps of ten from 0. 1 /urn to 1 .0 cm. 
A value of m = 3.5 is representative of interstellar grains (Mathis] 
|et al. 1977 1. A shallower slope of m = 2.5 is expected when 
grains coagulate to larger sizes ( [Natta et aL 2004||20()7] ), whereas 
a slope of m - 4.0 is expected when also fragmentation is taken 
into account ( Brauer et al. 2008 . (Dominik & Dullemond 2008] ). 
The different model parameters are summarised in Table |6] 

The resulting SEDs from six models, with Omax varying and 
the other parameters kept fixed, is shown in Fig.|6] In these mod- 
els, 7?out was fixed to 300 AU and the scale height was kept fixed 
at the same value in all models to show only the effect of varying 
fmax- Strong variations are seen in all wavelength regimes, from 
the near-infrared through the mm. At wavelengths A < 2 yum, 
grain size distributions without grains larger than 1 /im give such 
a high opacity that the central star is significantly reddened. In 
the mid- and far infrared, the flux drops with increasing maxi- 
mum grain size. The (sub)mm part of the SED does not change 
appreciably unless grains with sizes of ~ 100 fim or larger are in- 
cluded. After that, the (sub)mm slope becomes shallower quite 
rapidly with increasing flmax- This figure also demonstrates that 
care must be taken when estimating the disc mass from the 
(sub)mm emission alone: even when the dust composition is 
kept the same, assuming a different grain size distribution may 
already change the opacity at 1 mm by an order of magnitude, 
which will give an equally large uncertainty in the mass estimate 
from an observed flux at that wavelength. 
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Fig. 6. Spectral energy distributions (SEDs) for models of a 5 x 
10"^ Mq disc with a varying grain size distribution, seen under 
an inclination / = 45°. The minimum grain size is 0.003 //m 
for all models, and the maximum grain size varies from 0. 1 yum 
(black curve) to 1.0 cm in steps of ten (dark blue, light blue, 
green, yellow, red). The inset shows a blow-up of the 10-/im 
features. Note that the SEDs for Omax = 0.1 ;um, 1.0 ytm, and 
10 yum completely overlap in the mm part. The sharp drop in flux 
for flmax = 100 /vm around 2000 yum is due to a drop in opacity 
at about /I = 2 x tt x fl,nax- 



4.2. 10-fim feature vs mm slope 

In Fig. [7] we plot the strength of the 10-fj.m feature vs the mm 
slope for different models. The strength of the 10-yum feature 
(Fio-/^cont)/^cont IS defined as in Furlan et al. (2006) and the mm 
slope a is determined between 1.0 and 3.0 mm. The main aim of 
this figure is to show the variation of the 10-//m-feature strength 
and mm slope with various parameters. While the quantitative 
details will depend on the specific dust and disc parameters used, 
the qualitative trends found in these figures should be robust. 

In each of the panels, the results for different maximum grain 
sizes are shown. The size of the triangles is an indication for the 
maximum grain size under consideration. A general trend is ob- 
served, in the sense that the models with only small grains end up 
in the upper right corner of the micron-vs-mm diagram (strong 
lO-yLfm feature and steep mm slope), the models which include 
grains of mm sizes or larger end up more to the lower left of 
the diagram (weak 10-fj.m feature and shallower mm slope), and 
those with grain sizes of up to 10 or 100 /im end up towards the 
upper left comer of the diagram (weak 10-yum feature and steep 
mm slope). A possible evolutionary sequence, in which the max- 
imum grain size in the disc gradually increases, is indicated by 
the arrows: first, the 10-yum feature becomes weaker and later, 
the mm slope becomes shallower. A test to check whether ra- 
dial variation of a^ax - larger grains closer to the star, where the 
densities are higher - did not show any significant difference. 

The models show the effect of the temperature and luminos- 
ity of the central star, on the strength of the 10-yum feature and 
the steepness of the mm slope. The left column shows the re- 
sults for a central star with Teff = 3000 K and L = 1 L© and 
the right column for Teff = 4000 K and L = 6 Lq. In Figs. 
and b, the power-law slope of the grain size distribution is var- 
ied from m = 2.5 to 3.0, 3.5, and 4.0. It appears that only grain 
size distributions with m = 2.5 produce completely flat 10-yum 



silicate features as well as mm slopes with a < 2.0, whereas 
grain size distributions with m - 4.0 never produce a mm slope 
with ff < 3.0. Furthermore, the strongest 10-jum features are only 
obtained with a central star of 4000 K and L - 6Lq. 

In Figs.|7|; and d, the power-law slope of the grain size dis- 
tribution is fixed to m - 3.5. The disc radius 7?out is varied 
between 100, 200, and 300 AU. This has a small effect on the 
strength of the 10-//m feature, particularly for Teff = 3000 K and 
«max = 0.1 jum. This can be understood in the sense that for a 
smaller disc with the same dust mass, a larger amount of mass is 
closer to the star and thus radiates in the infrared. The mm slope 
of the SED is practically unaffected. 

Figs. [7^ and f show the results for models in which the 
power-law slope of the grain size distribution was fixed to m = 
3.5, the disc outer radius to /?out = 300 AU, and for which the 
inclination / under which the disc is observed is varied. In most 
cases, the inclination has a limited effect on both the strength 
of the 10-yum feature and the mm slope of the SED. Only under 
very high inclination (e.g., 75°, where 90° is edge-on) does the 
10-yum feature appear in absorption (not shown). A similar effect 
is seen if the discs are more flaring than found in vertical hy- 
drostatic equilibrium: the 10-/im feature is primarily weakened, 
because of the enhanced extinction under most inclinations. 

A second set of models is run to investigate the effects of 
dust settling, i.e., the process in which larger grains fall to the 
disc midplane under the influence of gravity, while the smaller 
grains stay suspended in the disc atmosphere. As mentioned be- 
fore, Dullemond & Dominik ( 2008j ) found that a bimodal grain 
size distribution can explain variations in the strength of the 10- 
fim feature, but only under specific circumstances. They looked 
at grains that are mainly responsible for the 10-jum feature, in 
particular grains of 3 fim and of <s; 1 yum. To study the effect of 
the settling of larger grains, we ran a number of models with up 
to six different grain size distributions: grains with sizes between 
0.003 and 0.1 jum, between 0.1 and 1 yum, 1 and 10, 10 and 100, 
100 and 1000, and finally between 1000 and 10,000 fj.m. The de- 
gree of settling is given by a parameter s, varying between 0.25, 
0.50, 0.75, and 1.00, and is chosen to be different for each of 
the grain size distributions: the larger the grains, the larger the 
degree of settling. For example, if H denotes the self-consistent 
scale height, a settling parameter s = 0.75 indicates that: 

grains between 0.003 and 0. 1 jum are at H; 

grains between 0.1 and 1 yum are at 0.75 x H; 

grains between 1 and 10 fim are at 0.75^ x H; 

grains between 10 and 100 yum are at 0.75^ x H; 

grams between 100 and 1000 //m are at 0.75^* x H; 

grains between 1000 and 10,000 jum are at 0.75^ x H. 

Hence, a larger number for s indicates a smaller degree of 
settling and i = 1 .00 corresponds to no settling (all grains are at 
the self-consistent scale height). 

These models are run using the radiative transfer code 
MCMax (Min et al. 2009). MCMax and RADMC were bench- 
marked against the results of Pascucci et al.| ( 2004| and the dif- 
ferences in the resulting SEDs are minimal, with in particular 
the 10-yum features being practically indistinguishable (see the 
Appendix in Min et al. 2009 1. The results for the settling are 
shown in Fig.|8] Note that the strength of the 10-yum feature in- 
creases when settling is switched on, but does not increase sig- 
nificantly more when the value of s is decreased more (i.e., when 
the degree of settling is increased). The slope in the mm part of 
the SED is practically unaffected by the degree of settling. 

It can be concluded that a variation of the maximum grain 
size fltnax affects both the strength of the 10-yum feature and the 
steepness of the mm slope a: a larger maximum grain size yields 
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Fig. 8. Strength of the lO-j-im feature as a function of the settling 
parameter s and the maximum grain size flmax- The size of the 
points indicates Omax, which varies from 1 to 10, 100, 1000, and 
10,000 /jm. Note that a small value for s corresponds to a large 
degree of settling; s = 1.00 corresponds to no settling (all dust 
grains at the self-consistent scale height). All other disc param- 
eters are as for the standard model in Table |6] 

a weaker 10-;um feature and a shallower mm slope. This effect is 
robust against variations in the degree of settling s, which only 
affects the strength of the 10-//m feature. Only varying the max- 
imum grain size is, however, not enough to explain the range in 
10-/im features and mm slopes that is observed. Variations in the 
power-law index of the grain size distribution and the tempera- 
ture and luminosity of the central source are required as well to 
cover the full observed range. For example, only models with a 
central-star temperature of at least 4000 K yield a 10-/im feature 
with a strength of (Fio - Fcont)/^^cont ~ 1-0; models with rela- 
tively flat or shallow grain size distributions are required to get 
mm slopes of a < 2.0 (see also [Natta et al.|[2007 Ricci et al. 
|20T0l l. 

5. Discussion 

A correlation between the strength of the 10-/im silicate feature 
and the slope between 1 and 3 mm is observed in a sample of 
in total 31 T-Tauri stars. This seems to imply that, while in the 
upper layers of the hot inner disc small grains are still coagulat- 
ing to a few microns in size, in the cold mid-plane of the outer 
disc large grains of at least several millimetres in size are akeady 
being formed. One should keep in mind, though, that the 10-yum 
silicate feature only probes the hot surface layers of the inner 
disc. It is therefore possible, and even likely, that large grains 
or boulders exist in the mid-plane of the inner disc. The tenta- 
tive correlation seen in our observations show, however, that the 
inner and outer discs of young stellar objects do not evolve inde- 
pendently from each other. Furthermore, the gradual decline of 
the 10-yL(m feature as the mm slope becomes shallower implies 
that micron-sized grains remain present in the disc atmosphere 
as cm-sized pebbles are already forming in the mid-plane. 

The radiative-transfer programmes RADMC and MCMax 
were used to run a number of models with varying disc param- 
eters, stellar parameters, and power-law grain size distributions. 
The only parameter that affects the strength of the 10-//m fea- 
ture and the mm slope of the SED simultaneously is the max- 
imum grain size a^ax- a larger Omax gives both a weaker 10- 



l^m feature and a shallower mm slope. This result is very ro- 
bust against variations of the other parameters. There are, how- 
ever, a few striking results when the models are compared to 
the observations. Firstly, the strongest lO-fim silicate features, 
with {Fio - f'cont)/^^cont ~ 1-0, Can only be reproduced for a 
star with T^s - 4000 K and L - 6 Lq, whereas a tempearature 
of 3000 K and a luminosity of 1 Lq seem to be more common 



for observed T-Tauri stars (see Evans et al.^200 9 ). This may in 
part be a selection effect: the sources with the hottest central 
stars are the most likely to have strong mm fluxes, making it 
more likely for them to provide a good mm slope. A subsam- 
ple of the sources of Evans et al. (2009 ) which had enough data 
points in the SED to provide a decent model fit also gave on aver- 
age a higher T^ff than would be expected from the near-infrared 
colours alone (B. Merm, priv. comm.). The physical reason for 
the increase in strength of the lO-fim feature with temperature 
lies in the fact that the photosphere is still an important contin- 
uum source. When the temperature of the star is increased, the 
relative contribution from the photosphere becomes smaller and 
the peak/continuum ratio of the feature goes up. 

Secondly, the sources with (Fu) - Fcont)/^^cont ~ 0.0 as well 
as those with a mm slope a < 2.0 can only be reproduced if 
the power-law slope of the grain size distribution is as low as 
m = 2.5. A slope m = 4.0 does not allow the mm slope to get 
below 3.0. The strength of the 10-jum feature remains at least as 
large as 1.4, and even 1.6 in the case Teff = 4000 K. This ap- 
pears to indicate that, as grains are growing to mm and cm sizes, 
the power-law slope of the grain size distribution becomes shal- 
lower. It is, however, more likely that a power-law grain size 
distribution with a mimimum grain size of Gmin = 0.003 /urn 
no longer applies. Either the effective minimum grain size in- 
creases, as suggested by the work of Olofsson et al.| ( 2009[ l, or 
a different grain size distribution is more applicable, such as 
naturally obtained from models including fragmentation (e.g., 
Dul lemond & Donunik|2005] l. 

It is important at this point to note once more that the 
10-;um feature originates entirely from the disc atmosphere; 
(sub)micron-sized dust grains that are closer to the midplane 
are not observable due to the optical thickness of the infrared 
emission. The emission at (sub)mm wavelengths, on the other 
hand, is predominantly optically thin; hence, it may originate 
from the whole disc. In fact, the bulk of the (sub)mm emis- 
sion comes from the midplane, since the larger grains will set- 
tle to the midplane. To study the effect of the settling of larger 
(mm and cm-sized) grains, we ran a number of models with 
up to six different grain size distributions. The smallest grains 
(0.003 < fl < 0. 1 fj.m) followed the gas, while the increasingly 
larger grains were located increasingly closer to the midplane. It 
turns out that the strength of the IQ-fim feature always increases 
when settling is switched on, but does not increase significantly 
more when the value of the settling parameter s is decreased 
more (i.e., when the degree of settling is increased). This can 
be understood from the fact that a small degree of settling al- 
ready removes the larger grains from the disc atmosphere. Since 
the lO-yKm feature originates in the disc atmosphere, it does not 
matter where the larger grains reside, as long as they are not too 
close to the surface. The slope in the mm part of the SED is 
practically unaffected by the degree of settling. This is because 
the disc is optically thin at these long wavelengths and hence 
it does not matter where in the disc the large grains, which de- 
termine the brightness at long wavelengths, are located. While 
the absolute strength and slope depend on the detailed dust and 
disc model parameters, the trends found here should be robust. 



14 



Dave J. P. Lommen et al.: Correlating yum and mm observations in YSOs 



Further modelling and exploration of the large parameter space 
are left to a future paper. 

The observations indicate a clustering of the sources per 
star-forming region. The T-Tauri stars in Chamaeleon show the 
strongest 10-/im features, as well as the steepest mm slopes. The 
T-Tauri stars located in the Lupus 1 and Lupus 2 clouds have 
somewhat weaker 10-fj.m features and mm slopes that are a bit 
shallower than the Chamaelon sources. The sources located in 
the Taurus-Auriga star-forming region, finally, have the weak- 
est 10-//m features and the shallowest mm slope. Although the 
number statistics are too low to draw strong conclusions from 
this grouping per star-forming region, it is intersesting to hy- 
pothesise as to what might be the origin of this grouping. It 
could in principle be due to a selection effect. To check for this, 
the strengths of the 10-//m features of the eleven sources from 
Furlan et al.| (|2006[) included in this work were compared with 
those of the 72 Taurus sources not included in this work and 
the two samples were found to be statistically indistinguishable. 
Furthermore, after this paper was submitted, Ricci et al. ( 2010| l 
published mm slopes for a total of 21 sources in the Taurus- 
Auriga star-forming region. |Ricci et al.| found the slopes between 
1 and 3 mm all to be a < 2.7; the corresponding strengths of the 
IQ-fim features are (Fio - Fcont)/^cont < 0.75 for 17 of those 21 
sources ( Furlan et al.||2006) l. Thus, these sources mainly popu- 
late the lower left corner of the yum-vs-mm diagram. It is then 
tempting to attribute the clustering to an evolutionary sequence, 
with the Chamaeleon sources being the least evolved and the 
Taurus-Auriga sources the most evolved. If the evolution pro- 
ceeds equally fast in each star-forming region, Taurus-Auriga 
would then be the oldest star-forming region and Chamaeleon 
the youngest. However, ages of pre-main-sequence stars are no- 
toriously difficult to determine and there is a large spread of stel- 
lar ages within each region (see Table [TJ. Note that a chemical 
differentiation, with, e.g., Chamaeleon a region with more sili- 
cates than Taurus, would mainly have an effect on the strength of 
the 10-/im feature, whereas the mm slope of the SED is largely 
determined by the sizes of the grains in the disc alone ( Draine 
2006[ l. Also, sources are likely to go through their evolutionary 
stages at different rates, which may be determined by other pa- 
rameters such as initial conditions of the core. 

Our final sample only contained two cold discs for which the 
strength of the 10-/zm feature and the slope of the SED between 
1 and 3 mm could be obtained. These are T Cha, an evolved disc 
with no silicates in the spectrum that is not included in the anal- 
ysis, and CS Cha, a circumbinary disc. Several processes have 
been proposed which may cause the deficiency of hot dust in the 
inner disc, such as photo-evaporation, the clearing out of the in- 
ner disc by a planet, or grain growth up to mm sizes and larger 
If grain growth were the dominating reason for the decrease in 
infrared flux, one would expect the cold disc to show up in the 
lower part of the 10-//m-feature vs mm-slope diagram. However, 
CS Cha is located rather in the centre of the diagram and it would 
appear that grain growth is not the main reason for the lack of 
infrared emission. Indeed, since CS Cha harbours a binary, it 
may be the case that the inner disc is cleared out due to binary 
interactions, with a "normal" outer disc remaining. 

6. Conclusions 

We observed five binary systems and 35 single T-Tauri stars (of 
which one turned out to harbour two sources as well) at infrared 
and mm wavelengths. New Spitzer IRS spectra of the wavelength 
region containing the 10-yum silicate feature are presented for 13 
sources and the slopes in the mm regime of the SED are ob- 



tained, also for 13 sources. Furthermore, dust disc masses are 
obtained for 21 new sources, including strict 3cr upper limits for 
nine of the sources. Combining our new observations with data 
from the literature, a tentative correlation between the strength 
of the 10-yum feature and the mm slope is found for a sample of 
in total 31 T-Tauri stars. This seems to suggest that the inner and 
outer discs do not evolve independently. 

Modelling protoplanetary discs with varying geometries and 
grain size distributions indicates that grain growth alone cannot 
explain the observed evolution in the strength of the lO-yum fea- 
ture and the mm slope of the SED. It would rather seem that as 
the maximum grain size increases, also the power-law slope of 
the grain size distribution becomes shallower. Ricci et al. ( 2010| l 
also found that for the majority of the sources in their sample a 
grain size distribution flatter than that of the interstellar medium 
is required (see also|Natta et al.|2007|). This is an indication that 
a power-law grain size distribution with a fixed minimum grain 
size is too simple a representation of the dust in protoplanetary 
discs. 

The sample contains only one cold or transitional disc for 
which the mm slope in the SED and the strength of the 10-jum 
silicate feature could be obtained properly. This source, CS Cha, 
ends up rather in the centre of the IQ-fim vs mm-slope diagram, 
indicating that grain growth is probably not the source of the 
removal of dust from the inner disc. Transitional discs are often 
relatively weak in the (sub)mm regime and a new generation of 
telescopes is required to determine the mm slope for more of 
these sources. 

Although all 10-;um features used in this work were obtained 
using the Spitzer Space Telescope, the mm fluxes and hence 
the mm slopes were determined with a large number of differ- 
ent telescopes. As this potentially leads to systematic effects, it 
will be necessary to observe a statistically large enough sam- 
ple with only one or maybe two (sub)mm telescopes. Ideally, 
about a dozen sources per star-forming regions for at least four 
regions should be observed with the same telescope. However, 
due to the relative weakness of the sources under considera- 
tion, the scarcity of available telescope time, and the spread of 
young star-forming regions over the sky, this is currently hard 
to achieve. The Atacama Large Millimeter/submillimeter Array 
(ALMA) will, upon completion, be able to study an order of 
magnitude more sources than presented here, providing us with 
the statistics to study interesting relations such as the grouping of 
YSOs by star-forming region in the jum-vs-mm diagram. ALMA 
will also be able to resolve most of the sources studied here. 

This work has mainly focused on observations at 1 and 
3 mm, at which grains with sizes of the order of 1 cm can be stud- 
ied (Draine 2006). The detection of larger grains, or pebbles, in 
protoplanetary discs, requires observations at cm wavelengths. 
However, there may be other sources of emission at cm wave- 
lengths, such as magnetic fields or stellar winds, and it is nec- 
essary to monitor sources over extended periods of time to dis- 
entangle the different emission mechanisms ( [Wilner et al.|2005j 
[Lommen et"aL]|2009l l. The maximum bandwidth of the ATCA 
was recently improved by a factor of 16 with the implementa- 
tion of the Compact Array Broadband Backend (CABB). This 
will for the first time allow the survey of large samples of young 
stellar objects in the southern hemisphere at cm wavelengths. 

Thus, in five to ten years from now we will not only be study- 
ing how grain growth occurs in protoplanetary discs, but we will 
actually be able to pinpoint the locations in the discs where the 
formation of pebbles and boulders, the precursors to planetes- 
imals and planets, is taking place. Furthermore, on-going ob- 
servations at cm wavelengths will give us a stronger handle on 
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the processes taking place in and around the discs, affecting the 
growth of grains and the formation of planets. 
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Appendix A: Observations 

In Lupus, 15 single sources and the binary IK Lup & Sz 66 
were observed with the Submillimeter Array (SMA) at about 
1 mm. Furthermore, the binaries VV CrA and S CrA and the 
single source DG CrA in Corona Australis were observed with 
the SMA at about 1 mm. A full log of the SMA observations is 



given in Table A. 1 The results of these observations are shown 



in Table B.l and Fig. B.l 



With the Australia Telescope Compact Array (ATCA), 15 
sources spread over the constellations Lupus, Vela, Corona 
Australis, and Chamaeleon were observed at 3 and 7 mm. A log 



of the ATCA observations is given in Table A. 2 

The Combined Array for Research in Millimeter-wave 
Astronomy (CARMA) was used to observe eleven single sources 
and one binary in Serpens at 1 and 3 mm. The log of these ob- 
servations is presented in Table A. 3 the full results are given in 
Table iJl 

Finally, eight sources in Serpens were observed with the 
(Very Large Array) VLA at 7 mm and at 1.3, 3.6, and 6 cm. 



A full log of those observations is given in Table A. 4 



Appendix B: Resuits 

The complete results of the SMA observations are shown in 
Table |B.l| The amplitude as a function of (m, v) distance is plot- 



ted in Fig. B.l 



The complete results of the ATCA observations are shown 
in Table IB. 21 Several sources were observed at the same wave- 
length more than once. The data for those sources were co-added 
in the (m, v) plane to improve the signal-to-noise ratio. The re- 
sulting fluxes or 3cr upper limits are presented in Table |3] in the 
main text. 

The binary VV CrA was not resolved with the ATCA at 3 or 
7 mm using natural weighting, which is optimised for sensitivity. 
However, using uniform weighting, which is optimised for reso- 
lution, the binary could be resolved at 3 mm. The map is shown 



in Fig. B.2 The amplitude as a function of (m, v) distance of the 



B.3 



B.4 



sources detected with the ATCA at 3 mm is plotted in Fig 
that of the sources detected with the ATCA at 7 mm in Fig 

The complete results of the CARMA observations are shown 
in Table |B.3| Data from tracks that could not be properly cal- 
ibrated due to a weak gain calibrator are not included. Some 
sources were observed twice, once in the C and once in the D 
configuration. If good data were obtained in both occasions, they 
were co-added in the (m, v) plane to improve the signal-to-noise 
ratio. The resulting fluxes or 3cr upper limits are presented in 
Table [3] in the main text. The amplitude as a function of (m, v) 
distance of the sources detected with CARMA at 1 mm is plot- 
ted in Fig. B.5 that of the sources detected with CARMA at 
3 mm in Fig.|B.6| 
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Fig. B.2. Image of VV CrA, observed at 3.2 mm on 2 and 3 
August 2008. The offsets are with respect to the phase centre, 
which is located at 19:03:06.7, -37:12:49.7. The contours are at 
2, 4, 6, ... times the rms of 2.0 mJy/bm; negative contours are 
dashed. 



The complete results of the VLA observations are shown in 
Table 
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Table A.l. Overview of the SMA observations. 



Obs. date 


Wavelengths 
(mm) 


Config. 


Target sources" 


Gain cahbrator(s) 


Flux cal. 


Notes 


20080419 


1.30, 1.36 


compact 


EX Lup, RX J1615.3-3255 
kY Lup, lolUiy.j /-392214. / 
161159.81-382338.5, Sz 111 
Sz91,Sz 96 


1517-243, 1626-298 


Mars 


PWV < 4.0 mm 


20081001 


1.31, 1.37 


compact 


VV CrA, S CrA, DG CrA 


1924-292 


CaUisto 


PWV < 2.5 mm 


20090507 


1.30, 1.36 


compact 


HM Lup, MY Lup, 
Sz 102, IKLup&Sz66 
Sz 73, Sz 74 
Sz 76, Sz 77 


1626-298 


Callisto 





^ hi the case of SSTc2d names, only the coordinates (in J2000) are shown. 
Table A.2. Overview of the ATCA observations. 



Obs. date 


Wavelengths 
(mm) 


Config. 


Target sources 


Gain calibrator(s) 


Flux cal. 


Notes 


20080712 


3.14,3.21 


H214 


HBC 556, HBC 557, HBC 559 


0745-330 


Mars 


Mostly nice weather 


20080713 


6.67, 6.99 


H214 


HBC 556, HBC 557, HBC 559 


0826-373 


Uranus 


Weather worsening 


20080728 


6.67, 6.99 


H214 


HBC 559 


0826-373 


Uranus 


Mostly nice weather 


20080729 


3.14, 3.21 


H214 


HBC 559 


0745-330 


Uranus 


Mostly nice weather 








SZ Cha, Sz 32 


1057-797 






20080730 


3.14,3.21 


H214 


SZ Cha 


1057-797 


Uranus 


Weather less than perfect 


20080801 


6.67, 6.99 


H214 


Sz 111, RY Lup 


1600-44 


Uranus 


Weather improving 








RX J1615.3-3255 


1622-297 












VV CrA, S CrA, DG CrA 


1921-293 






20080802 


3.14,3.21 


H214 


SZCha 


1057-797 


Uranus 


Weather improving 








Sz 111, RY Lup 


1600-44 












RX J1615.3-3255 


1622-297 












VV CrA, S CrA, DG CrA 


1933-400 






20080803 


3.14,3.21 


H214 


MY Lup 


1600-44 


Uranus 


Mostly nice weather 








VV CrA, S CrA, DG CrA 


1933-400 






20080804 


3.14, 3.21 


H214 


Sz 65 & Sz 66 


1622-297 


Uranus 


Weather less than perfect 


20080805 


6.48, 6.78 


H214 


MY Lup, IM Lup 


1600-44 


Uranus 


Weather improving 
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Wavelength 
(mm) 


Target source 


Conflg. 


Obs. date(s) 


Gain calibrator(s) 


Flux cal. 


Notes 


1.33 


EC 82 


C 


20080424 


1743-038 


3c454.3 


Gain calibrator too weak 




EC 90 


C 


20080424 


1743-038 


3c454.3 


Gain calibrator too weak 




SSTc2d J182900.88+002931.5 


c 


20080424, 20080426 


1743-038 


3c454.3 


Gain calibrator too weak 




IRAS 18268-0025 


c 


20080426 


1743-038 


3c454.3 


Gain calibrator too weak 




CoKu Ser-G3 


c 


20080426 


1743-038 


3c454.3 


Gain calibrator too weak 




SSTc2d J182858.08+001724.4 


c 


20080426 


1743-038 


3c454.3 


Gain calibrator too weak 




SSTc2d J182850.20+000949.7 


c 


20080426 


1743-038 


3c454.3 


Gain cahbrator too weak 




SSTc2d J182944.10+003356.1 


c 


20080426 


1743-038 


3c454.3 


Gain cahbrator too weak 




VVSer 


c 


20080426 


1743-038 


3c454.3 


Gain cahbrator too weak 




SSTc2d J182936.19+004216.7 


c 


20080501 


1743-038 


3c454.3 


Gain cahbrator too weak 




GSC 00446-00153 


c 


20080501, 20080518 


1743-038, 1751-1-096 


3c454.3 




1.33 


EC 82 


D 


20080620 


175U096 


3c273 






EC 90 


D 


20080620 


175U096 


3c273 






SSTc2d 1182900.88-^002931.5 


D 


20080620 


175U096 


3c273 






IRAS 18268-0025 


D 


20080620 


175U096 


3c273 






CoKu Ser-G3 


D 


20080620 


1751-1-096 


3c273 






VVSer 


D 


20080620 


1751-1-096 


3c273 






SSTc2d J182858.08-H001724.4 


D 


20080620 


1751-1-096 


3c273 






SSTc2d Jl 82850.20-1-000949.7 


D 


20080620, 20080704 


1751-1-096 


3c273 






EC 97 


D 


20080704 


1751-^096 


3c279 






SSTc2d 1182944.10-^003356.1 


D 


20080704 


175U096 


3c279 






SSTc2d 1 1 82936. 1 9-^0042 1 6.7 


D 


20080704 


1751-1-096 


3c279 






GSC 00446-00153 


D 


20080704 


1751-1-096 


3c279 




3.15 


EC 82 


C 


20080521 


1743-038 


3c273 


Gain calibrator too weak 




EC 90 


C 


20080521 


1743-038 


3c273 


Gain cahbrator too weak 




SSTc2d J182900.88-H002931.5 


C 


20080521 


1743-038 


3c273 


Gain cahbrator too weak 




IRAS 18268-0025 


C 


20080521 


1743-038 


3c273 


Gain cahbrator too weak 




SSTc2d Jl 82858.08-1-001724.4 


C 


20080521 


1743-038 


3c273 


Gain calibrator too weak 




SSTc2d J182850.20-I-000949.7 


C 


20080521 


1743-038 


3c273 


Gain calibrator too weak 




VV Ser 


C 


20080521 


1743-038 


3c273 


Gain calibrator too weak 


3.15 


EC 82 


D 


20080618 


175U096 


3c279 






EC 90 


D 


20080618 


175U096 


3c279 






SSTc2d 1182900.88-^002931.5 


D 


20080618 


175U096 


3c279 






IRAS 18268-0025 


D 


20080618 


1751-1-096 


3c279 






CoKu Ser-G3 


D 


20080618 


1751-1-096 


3c279 






VVSer 


D 


20080618 


1751-1-096 


3c279 






SSTc2d J182858.08-H001724.4 


D 


20080619 


1751-1-096 


3c273 






SSTc2d 1182850.20-F000949.7 


D 


20080619 


1751-^096 


3c273 






SSTc2d 1182944.10-^003356.1 


D 


20080619 


1751-F096 


3c273 






SSTc2d 1 1 82936. 1 9-^0042 1 6.7 


D 


20080619, 20080622 


1751-1-096 


3c273 






GSC 00446-00153 


D 


20080622 


1751-1-096 


3c273 





Table A.4. Overview of the VLA observations. 



Obs. date 


Wavelengths 
(mm) 


Config. 


Target sources" 


Gain calibrator(s) 


Flux cal. 


Notes 


20080310 


6.92, 6.93 


C 


CoKu Ser G3, EC 82 
182900.88-1-002931.5, EC 90, VV Ser, 
EC 97, 182850.20-1-000949.7 


1824-1-013 


3c286 


Clouds forming 


20080311 


13.4 


C 


CoKu Ser G3, 182900.88-1-002931.5, 
EC 90, VV Ser, EC 97 


1851-1-005 


3c286 


Clear sky 


20080313 


13.4 


C 


182850.20-^000949.7, 
182909.80-^003445.9, EC 82 


1851-1-005 


3c286 


Clouds forming 


20080314 


35.5 


C 


EC 82, EC 90, 182950.20-1-000949.7, 

EC 97, VV Ser, 
CoKu Ser G3, 182900.88-1-002931.5, 


1804-1-010 


3c286 


Clear sky; high winds 


20080315 


61.4, 62.0 


C 


EC 82, EC 90, 182950.20-1-000949.7, 


1804-1-010 


3c286 


Clear sky; high winds 



EC 97, 

CoKu Ser G3, 18290088-1-0029315, 
VV Ser, 18285020-1-0009497 



" In the case of SSTc2d names, only the coordinates (in J2000) are shown. 
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Table B.l. Complete results of SMA observations at 1.3 nun. 



Obs. date 


Effective 


Target source" 


Continuum flux" 


rms*^ 


Gaussian size 


RA" 


Dec" 




wavelength 




(P) (G) 




(arc sec) 


(J2000) 


(J2000) 




(mm) 




(mJy) 


(mjy/bm) 








20080419 


1.34 


EXLup 


19.3 21.3 


4.0 


1.00 ± 1.34 


16 03 05.0 


-40 18 20.1 






RXJ1615. 3-3255 


131.8 169.1 


3.9 


1.53 ±0.13 


16 15 20.2 


-32 55 05.3 






RYLup 


78.3 89.0 


5.0 


1.14 ±0.30 


15 59 28.4 


-40 21 51.4 






161029.57-392214.7 


< 13. l"" 


4.4 


- 


16 10 29.6 


-39 22 14.4 






161159.81-382338.5 


< 12.7" 


4.2 


- 


16 11 59.8 


-38 23 38.0 






Sz 111 


49.3 52.5 


4.2 


0.78 ± 0.67 


16 08 54.7 


-39 37 43.6 






Sz91 


< 13.0" 


4.3 


- 


16 07 11.6 


-39 03 47.1 






Sz 96 


< 12.6" 


4.2 




16 08 12.6 


-39 08 33.3 


20081001 


1.35 


VVCrA 


349.6 367.0 


5.1 


0.9(1 ± 0.06 


19 03 06.8 


-37 12 49.3 






SCrA 


301.4 322.2 


3.5 


0.91 ± 0.07 


19 01 08.6 


-36 57 20.6 






DGCrA 


< 6.6" 


2.2 


- 


19 01 55.2 


-37 23 40.5 


20090507 


1.34 


IKLup 


29.4 29.4 


2.8 


(unresolved) 


15 39 27.8 


-34 46 17.8 






Sz 66 


< 8.3" 


2.8 


- 


15 39 28.3 


-34 46 18.0 






HMLup 


< 10.2" 


3.4 




15 47 50.6 


-35 28 35.3 






Sz73 a 


16.2 


2.9 




15 47 57.0 


-35 24 35.9 






Sz 73 b 


15.8 


2.9 




15 47 57.1 


-35 14 40.0 






Sz 74 


15.1 15.1 


3.0 


(unresolved) 


15 48 05.3 


-35 15 53.8 






Sz76 


12.4 12.4 


3.3 


(unresolved) 


15 49 30.8 


-35 49 51.2 






Sz77 


< 9.5" 


3.2 




15 51 47.0 


-35 56 44.1 






MYLup 


56.4 66.1 


3.4 


1.43 ±0.51 


16 00 44.5 


-41 55 31.5 






Sz 102 


< 11.4" 


3.8 




16 08 29.7 


-39 03 11.0 



^ In the case of SSTc2d names, only the coordinates (in J2000) are shown. 
Continuum flux and position are from fits in the (u, v) plane. For sources that were detected at 3cr, both the point-source flux (P) and the 
integrated flux for a Gaussian (G) are shown. For sources that were not detected, the coordinates of the phase centre are quoted. 
Calculated from the cleaned image. 

" Quoted value is 3cr upper limit. 

^ No circular Gaussian could be fit to the source in the {u, v) plane. 
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Obs. date 


Effective 


Target source 


Continuum flux^ 


rms" 


Gaussian size 


RA^ 


Dec" 




wavelength 




(P) 


(G) 




(arcsec) 


(J2000) 


(J2000) 




(mm) 




(mJy) 




(mjy/bm) 








20080712 


3.18 


HBC 556 


< 3.7'^ 




1.2 


- 


8 10 31.3 


-36 01 46.5 






HBC 557 


< 3.2'^ 




1.1 


- 


8 12 47.4 


-36 19 18.0 






HBC 559 


< 3.3'^ 




1.1 


- 


8 13 56.4 


-36 08 02.1 


20080713 


6.83 


HBC 556 


<0.7'= 




0.2 


- 


8 10 31.6 


-36 01 46.5 






HBC 557 


<0.4'^ 




0.1 


- 


8 12 47.7 


-36 19 18.0 






HBC 559 


<0.6'^ 




0.2 


- 


8 13 56.8 


-36 08 02.1 


20080728 


6.85 


HBC 559 


<0.3'^ 




0.1 


- 


8 13 56.8 


-36 08 02.1 


20080729 


3.18 


HBC 559 


< 4.3'= 




1.4 


- 


8 13 51.0 


-36 08 02.1 






SZ Cha 


< 3.0*= 




1.0 


- 


10 58 10.0 


-77 17 17.6 






Sz32 


<2.9'= 




1.0 


- 


11 09 48.0 


-76 34 26.0 


20080730 


3.18 


SZ Cha 


3.4 


5.8 


0.5 


1.95 ±0.63 


10 58 16.6 


-77 17 17.0 


20080801 


6.85 


Sz 111 


< 0.6'= 




0.2 


- 


16 08 53.8 


-39 37 43.1 






RYLup 


<0.6'= 




0.2 


- 


15 59 27.5 


-40 21 51.2 






RXJ1615.3-3255 


<0.5'= 




0.2 


- 


16 15 19.4 


-32 55 05.0 






VVCrA 


3.4 


3.4 


0.2 


(unresolved) 


19 03 06.8 


-37 12 49.3 






SCrA 


3.7 


5.4 


0.2 


3.12 ±0.81 


19 01 08.6 


-36 57 20.3 






DG CrA 


< 0.6^ 




0.2 


- 


19 01 54.4 


-37 23 40.5 


20080802 


3.18 


SZ Cha 


< 2.9'= 




1.0 


- 


10 58 15.5 


-77 17 17.6 






Sz 111 


5.7 


e 


0.7 


- 


16 08 54.6 


-39 37 53.3 






RYLup 


< 2.3'= 




0.8 


- 


15 59 28.0 


-40 21 51.2 






RX J1615.3-3255 


6.8 


_d 


0.6 


- 


16 15 20.2 


-32 55 05.6 






VV CrA a 


31.0 


_e 


1.2 


_e 


19 03 06.8 


-37 12 49.8 






VVCrAb 


25.1 uniform 


_e 




_e 


19 03 06.9 


-37 12 48.3 






SCrA 


22.0 


_d 


1.1 


- 


19 01 08.6 


-36 57 20.2 






DG CrA 


< 3.0"= 




1.0 




19 01 54.9 


-37 23 40.5 


20080803 


3.18 


MYLup 


8.7 


_d 


0.4 




16 00 44.5 


-41 55 31.2 






VV CrA a 


23.0 


_e 


1.9 




19 03 06.8 


-37 12 49.9 






VV CrA b 


21.9 


_e 




■ 


19 03 06.9 


-37 12 48.4 






SCrA 


24.9 


_d 


1.9 




19 01 08.6 


-36 57 20.6 






DG CrA 


<4.3'= 




1.4 




19 01 54.9 


-37 23 40.5 


20080804 


3.18 


Sz 65 


3.4 


_d 


0.4 




15 39 27.7 


-34 46 17.6 






Sz 66 


2.2 


_d 






15 39 28.2 


-34 46 17.9 


20080805 


6.65 


MYLup 


L3 


2.9 


0.1 


4.75 ± 1.31 


16 00 44.6 


-41 55 31.5 






IMLup 


2.2 


2.2 


0.2 


(unresolved) 


15 56 09.2 


-37 56 06.0 



" Continuum flux and position are from fits in the (m, v) plane. For sources that were detected at 3cr, both the point-source flux (P) and the 

integrated flux for a Gaussian (G) are shown. For sources that were not detected, the coordinates of the phase centre are quoted. 

Calculated from the cleaned image. 
■= Quoted value is 3cr upper limit. 

No circular Gaussian could be fit to the source in the («, v) plane. 
^ The two components could not be separated with circular Gaussian fits in the («, v) plane. One circular Gaussian was fitted to the binary, 

yielding a flux of 69.5 mJy and a size of 2.56 ± 0.21 arcsec on 2 August 2008 and a flux of 44.2 mJy and a size of 2.28 ± 0.38 arcsec on 3 

August 2008. 
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Obs. date 


Effective 


Target source" 


Continuum flux" 


rms*^ 


Gaussian size 


RA" 


Dec" 




wavelength 




(P) (G) 




(arc sec) 


(J2000) 


(J2000) 




(mm) 




(mJy) 


(mjy/bm) 








20080518 


1.33 


GSC 00446-00153 


43.6 50.5 


6.0 


0.46 ± 0.32 


18 30 06.2 


-1-00 42 33.6 


20080618 


3.15 


EC 82 


< 2.9'' 


1.0 


- 


18 29 56.8 


-1-01 14 46.0 






EC 90 


11.3 11.6 


1.0 


1.69 ± 0.97 


18 29 57.8 


-1-01 14 06.9 






182900.88+002931.5 


3.4 3.9 


0.6 


(unresolved) 


18 29 00.9 


-1-00 29 31.7 






IRAS 18268-0025 


< 1.9'' 


0.6 


— 


18 29 28.2 


-00 22 57.1 






CoKu Ser G3 


< 1.6" 


0.5 


— 


18 29 01.8 


-1-00 29 54.6 






VV Ser 


< 1.8" 


0.6 


— 


18 28 47.9 


-1-00 08 40.0 


20080619 


3.15 


182858.08-1-001724.4 


< 1.9" 


0.6 


— 


18 28 58.1 


-1-00 17 24.4 






182850.20-1-000949.7 


< 1.9" 


0.6 


— 


18 28 50.2 


-1-00 09 49.7 






182944.10-1-003356.1 


< 1.7" 


0.6 


— 


18 29 44.1 


-1-00 33 56.1 






182936.19-1-004216.7 


< 2.7" 


0.9 


— 


18 29 36.2 


-1-00 42 16.7 


20080620 


1.33 


EC 82 


< 15.7'' 


5.2 


— 


18 29 56.8 


-1-01 14 46.0 






EC 90 


91.8 91.7 


10.4 


(unresolved) 


18 29 57.7 


-1-01 14 07.0 






182900.88-1-002931.5 


< 16.4" 


5.5 


— 


18 29 00.9 


-1-00 29 31.6 






IRAS 18268-0025 


< 15.6" 


5.2 


— 


18 29 28.2 


-00 22 57.1 






CoKu Ser G3 


< 17.1" 


5.7 


— 


18 29 01.8 


-1-00 29 54.6 






VV Ser 


< 14.8" 


4.9 


— 


18 29 47.9 


-1-00 08 40.0 






182858.08-1-001724.4 


< 24.6" 


8.2 


— 


18 28 58.1 


-1-00 17 24.4 






182850.20-1-000949.7 


< 98.2" 


32.7 




18 28 50.2 


+00 09 49.7 


20080622 


3.15 


182936.19-^004216.7 


< 5.3" 


1.8 




18 29 36.2 


-1-00 42 16.7 






GSC 00446-00153 


6.8 8.1 


1.0 


2.71 ± 1.58 


18 30 06.3 


-1-00 42 34.2 


20080704 


1.33 


EC 97 


< 23.3" 


7.8 




18 29 58.2 


-1-01 15 22.0 






182850.20-^000949.7 


< 22.8" 


7.6 




18 28 50.2 


-1-00 09 49.7 






182944.10-^003356.1 


< 15.0" 


5.0 




18 29 44.1 


-1-00 35 56.1 






182936.19-1-004216.7 


< 8.7" 


2.9 




18 29 36.2 


-1-00 42 16.7 






GSC 00446-00153 


90.8 97.6 


2.8 


0.67 ± 0.23 


18 30 06.2 


-1-00 42 33.6 



^ In the case of SSTc2d names, only the coordinates (in J2000) are shown. 
Continuum flux and position are from fits in the («, v) plane. For sources that were detected at 3cr, both the point-source flux (P) and the 
integrated flux for a Gaussian (G) are shown. For sources that were not detected, the coordinates of the phase centre are quoted. 
Calculated from the cleaned image. 

" Quoted value is 3cr upper limit. 
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Fig. B.5. Amplitude as a function of {u, v) distance for sources detected with CARMA at 1 mm. 
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Fig. B.6. Amplitude as a function of (m, v) distance for sources detected with CARMA at 3 mm. 
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Table B.4. Complete results of VLA observations at 6.9 mm and 1.3, 3.5, and 6.2 cm. 



Obs. date 


Effective 


Target source" 


Continuum flux" 


rms'^ 


Gaussian size 


RA" 


Dec" 




wavelength 




(Peak) (Integ.) 




(arc sec) 


(J2000) 


(J2000) 




(mm) 




(mJy) 


(mjy/bm) 








20080310 


6.93 


CoKu Ser G3 


< 1.2'' 


0.4 


— 


18 29 01.8 


+00 29 54.5 






182900.88+002931.5 


< 0.7'' 


0.2 


— 


18 29 00.9 


+00 29 31.5 






EC 90 


< l.O"* 


0.3 


— 


18 29 57.7 


+01 14 05.7 






VVSer 


< 0.7'' 


0.2 


— 


18 28 47.9 


+00 08 39.8 






EC 97 


< 0.6'' 


0.2 


— 


18 29 58.2 


+01 15 21.7 






182850.20+000949.7 


< 0.6'' 


0.2 


— 


18 28 50.2 


+00 09 49.6 






EC 82 


< 0.5'' 


0.2 


— 


18 29 56.9 


+01 14 46.4 


20080311 


13.4 


CoKu Ser G3 


8.66 9.81 


0.06 


1.17 X 0.91 


18 29 01.8 


+00 29 54.8 






182900.88+002931.5 


< O.IT^ 


0.06 


— 


18 29 00.9 


+00 29 31.5 






EC 90 


< 0.24'' 


0.08 


— 


18 29 57.7 


+01 14 05.7 






VVSer 


< 0.17'' 


0.06 


— 


18 28 47.9 


+00 08 39.8 






EC 97 


< 0.15'' 


0.05 


— 


18 29 58.2 


+01 15 21.7 


20080313 


13.4 


182850.20+000949.7 


< 0.20'' 


0.07 


— 


18 28 50.2 


+00 09 49.6 






182909.80+003445.9 


< 0.23'' 


0.08 


— 


18 29 09.8 


+00 34 45.8 






EC 82 


< 0.28'' 


0.09 


— 


18 29 56.9 


+01 14 46.4 


20080314 


35.5 


EC 82 


< 0.09'' 


0.03 


— 


18 29 56.0 


+01 14 49.9 






EC 90 


< 0.09'' 


0.03 


— 


18 29 56.0 


+01 14 49.0 






EC 97 


< 0.09" 


0.03 


— 


18 29 56.0 


+01 14 49.0 






CoKu Ser G3 


1.11 1.11 


0.03 


2.49 X 2.35 


18 29 01.8 


+00 29 54.7 






182900.88+002931.5 


< 0.10" 


0.03 


- 


18 29 05.0 


+00 29 44.0 






VV Ser 


0.14 0.17 


0.03 


3.25 X 2.31 


18 28 47.9 


+00 08 40.1 






182850.20+000949.7 


< 0.08'' 


0.03 




18 28 40.0 


+00 09 13.0 


20080315 


61.7 


EC 82 


< 0.09" 


0.03 




18 29 57.0 


+01 14 40.0 






EC 90 


< 0.09" 


0.03 




18 29 57.0 


+01 14 40.0 






EC 97 


< 0.09" 


0.03 




18 29 57.0 


+01 14 40.0 






CoKu Ser G3 


0.91 1.00 


0.04 


(unresolved) 


18 29 01.9 


+0029 54.9 






182900.88+002931.5 


< 0.12" 


0.04 




18 29 07.0 


+00 32 09.0 






VV Ser 


< 0.14" 


0.05 




18 29 47.0 


+00 09 11.0 






182850.20+000949.7 


< 0.14" 


0.05 




18 29 47.0 


+00 09 11.0 



" In the case of SSTc2d names, only the coordinates (in J2000) are shown. 

Continuum flux and position are obtained using the AIPS task JMFIT, which fits elliptical Gaussians to the cleaned image. For sources that 

were detected at 3cr, both the peak (Peak) and the integrated (Integ.) flux are shown. For sources that were not detected, the coordinates of the 

phase centre are quoted. 

Calculated from the cleaned image. 
" Quoted value is 3cr upper limit. 



